KEK-TH-648 

SNS-PH/99-16 

hep-ph/9912260 



Supersymmetry versus precision experiments 

revisited 

q ', Gi-Chol ChoQ 1 ^ 2 ) and Kaoru Hagiwara 1 ) 

o ; 

1 - ) Theory Group, KEK, Tsukuba, Ibaraki 305-0801, Japan 
2 ^Scuola Normale Superiore, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy 

a 

Abstract 

We study constraints on the supersymmetric standard model from the updated elec- 



O 



ON 



troweak precision measurements — the Z-pole experiments and the PP-boson mass mea- 



surements. The supersymmetric-particle contributions to the universal gauge-boson- 
propagator corrections are parametrized by the three oblique parameters Sz, T z and 



m w . The oblique corrections, the Zqq and Zll vertex corrections, and the vertex and box 
corrections to the /i-decay width are separately studied in detail. We first study individual 
contribution from the four sectors of the model, the squarks, the sleptons, the supersym- 
metric fermions (charginos and neutralinos) , and the supersymmetric Higgs bosons, to 
the universal oblique parameters, where the sum of individual contributions gives the 
total correction. We find that the light squarks or sleptons, whose masses just above the 
present direct search limits, always make the fit worse than that of the Standard Model 
(SM), whereas the light charginos and neutralinos generally make the fit slightly better. 
The contribution from the supersymmetric Higgs sector is found small. We then study 
the vertex/box corrections carefully when both the supersymmetric fermions (-inos) and 
the supersymmetric scalars (squarks and sleptons) are light, and find that no significant 
improvement over the SM fit is achieved. The best overall fit to the precision measure- 
ments are found when charginos of mass ~ 100 GeV with a dominant wino-component are 
present and the doublet squarks and sleptons are all much heavier. The improvement over 
the SM is marginal, however, where the total x 2 of the fit to the 22 data points decreases 
by about one unit, due mainly to a slightly better fit to the Z-boson total width. 
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1 Introduction 



The supersymmetric standard model has been the leading candidate for the the- 
ory beyond the Standard Model (SM) of elementary particles. The supersym- 
metry (SUSY), the symmetry between bosons and fermions, gives us an elegant 
solution |] to the hierarchy problem || , the stability of the electroweak scale 
(~ 100 GeV) against the more fundamental scale of physics, the Grand Unified 
Theory (GUT) scale (~ 10 16 GeV) where the three gauge interactions may be 
unified, or the Planck scale (~ 10 19 GeV) where the gauge interactions may be 
unified with the gravity interactions. It also offers us an attractive scenario that 
the electroweak symmetry breaking may occur radiatively Q if the top quark is 
sufficiently heavy but not too much (100 GeV <,m t <,200 GeV) 0. The obser- 
vation that the minimal supersymmetric standard model (MSSM) leads to the 
unification of the three gauge couplings at the 'right' scale, ~ 2 x 10 16 GeV, high 
enough for the proton longevity, if the supersymmetric particles exist with masses 
at or below the 1 TeV scale 0, and the discovery of the top quark [[F] in the 'right' 
mass range, jointly have made us take the MSSM as a serious candidate for the 
theory just beyond the present new-particle search front. 

Despite the high expectations, however, various efforts in search of an evi- 
dence for supersymmetric particles have so far been fruitless. There had been 
occasions when a short-lived experimental 'anomaly' such as the mismatch of a s 
determined from the low-energy experiments and that obtained from the Z-boson 
decays and/or the significantly large partial width T(Z — > bb) as compared to the 
SM prediction, were considered as a possible evidence for relatively light supersym- 
metric particles || |10|| . All such anomalies were short-lived and the precision 
electroweak data IJ] from the completed experiments at LEP1 and the Tevatron 
run-I, and the on-going experiments at SLC and LEP2 are all consistent with the 
SM predictions. Alongside, the direct search limits for the supersymmetric parti- 
cle masses have steadily risen with time. The only additional encouragement from 
the precision electroweak experiments is that the preferred range of the SM Higgs 



boson mass is consistent with the stringent upper bound jrj, [13], |Tj, [15|, |16| of the 
lightest Higgs boson mass in the supersymmetric theories. The good news is that 
the supersymmetric theories have not been ruled out yet, while the bad news is 
that studies of precision electroweak experiments give us only the lower bounds 
for the supersymmetric particle masses. 

There have recently been numerous attempts [[17], |18[ to identify constraints on 
the MSSM parameters from precision electroweak measurements. In this report, 
we present the first results of our comprehensive study of the constraints on the 
MSSM parameters from the electroweak precision measurements. 

We perform a systematic study of the MSSM parameter space by observing that 
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the concept of the universal gauge-boson-propagator corrections, or the oblique 
corrections [|1£], [21], |22|], is useful in the MSSM in the sense that they dom- 



inate the radiative corrections if either the supersymmetric scalars (squarks and 
sleptons) or the supersymmetric fermions (charginos and neutralinos) are all suffi- 
ciently heavy. In this limit, all the precision measurements on the Z-boson proper- 
ties can be parametrized by just two oblique parameters, Sz and T z , whereas the 
VT-boson mass itself makes the third oblique parameter. Here we closely follow 
the formalism developed in refs. [^, |25|, slightly modified to suit our MSSM 



studies. In this limit of the dominant oblique corrections, we can study in detail 
the contributions of the four sectors of the MSSM separately [p6|1 ; the squarks, the 
sleptons, the supersymmetric fermions (charginos and neutralinos), and the super- 
symmetric Higgs sector. We find that although the contributions from relatively 
light squarks and/or sleptons generally make the fit to the electroweak data worse 
than the SM fit, those from relatively light charginos and neutralinos can slightly 
improve the fit. 

In the second stage, we examine the case where both the supersymmetric scalars 
and fermions are light, by studying their contributions to the muon lifetime, A5 G , 
and to all the Z-boson decay amplitudes, Ag v L , Ag e Ll Ag R , Agl, Ag R , Agf, Ag R , 
Ag^J , Agl, Ag R , Ag b L , and Ag b R . Here Ag[, stands for the deviation of the non- 
oblique correction to the Zf a f a vertex from the reference SM prediction. Within 
our approximation of the MSSM Lagrangian, the amplitudes for the first two 
generation quarks and leptons are identical, and only the above 12 distinctive 
corrections appear. Because we find an indication that the existence of relatively 
light (~ 100 GeV) charginos and neutralinos can somewhat improve the SM fit to 
the electroweak data, we study the consequences of having light sleptons and/or 
squarks in addition to the light charginos and neutralinos. We generally find 
that the goodness of the fit worsens in such cases. Studies of the consequences 
of more specific models of the supersymmetry breaking mechanism, such as the 
supergravity models and the gauge-mediated supersymmetry breaking models, will 



be presented elsewhere [27 



The paper is organized as follows. In section 2, we present the formalism 
which we calculate all the electroweak observables in the MSSM. The dependences 
of those observables in terms of the SM parameters, m t , rriH SM , oc s = cn s (m z )y^, 



and a(m z ) are given in a compactly parametrized form [gj, gjj. The MSSM La- 



grangian and constraints on its parameters are summarized in section 3, and all 
the MSSM contributions to the one-loop functions are presented in the appendices. 
In section 4, we study oblique corrections from the squarks and sleptons (sec. 4.1), 
the supersymmetric Higgs bosons (sec. 4.2), and the supersymmetric fermions 
(charginos and neutralinos) (sec. 4.3), separately. In section 5, we study conse- 
quences of vertex and box corrections in addition to the gauge-boson-propagator 
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corrections: the vertex corrections due to the MSSM Higgs bosons (sec. 5.1), those 
from squarks and gluinos (sec. 5.2), those from squarks and charginos/neutralinos 
(sec. 5.3), and the vertex/box corrections from sleptons and charginos/neutralinos 
(sec. 5.4). In section 6, we examine quantitative significance of the MSSM con- 
tributions to the electroweak radiative corrections. Section 7 gives conclusions. 
Appendix A summarizes our notation |2S|| for the MSSM mixing matrices and the 



couplings among the mass eigenstates. Appendix B gives the MSSM contribution 
to the gauge-boson-propagator corrections and the oblique parameters S, T, U and 
R. Appendix C gives the MSSM contribution to the Z-boson decay amplitudes. 
Appendix D gives the MSSM contribution to the muon-decay amplitude. 

2 Electroweak observables in the MSSM 

In this section, we give all the theoretical predictions for the electroweak observ- 
ables which are used in our analysis. The experimental data of the Z-pole experi- 
ments and the J^-boson mass measurement are summarized in Table [|. Because 
the MSSM predicts slightly different Z-decay amplitudes into the r-leptons as 
compared to those into e or /x, we will use in the following analysis the data set 
which do not assume the e-/i-r universality. By removing the two entries, Re and 
y4.pg tj licit assume the lepton universality, there are 19 data points in the Table. 
The total \ 2 is obtained by using the correlation matrices of ref. JTT[ . 



The reference SM predictions are given in the table for an arbitrarily chosen 
set of the four inputs, m t = 175 GeV, m HsM = 100 GeV, l/a(m|) = 128. 90f] and 
a s = cn s (m z )y[s = 0.118. In the following, we present simple parametrizations 
of the SM predictions for the parameter sets around the above reference point in 
terms of the normalized variables; 

= m ,(GeV) - 175 

* 10 K J 

m HsM (GeV) 

h 100 K J 

l/aiml) - 128.90 

*° = "^09 ' (2 ' 1C) 

= q a (m z )Ms- Q.118 

0.003 v ; 

The SM predictions for an arbitrary set of the four input parameters (m t , rriH SM , 

a(m|), a s (m z )) are then obtained easily for all the electroweak observables. The 

* a(m 2 z ) is the running QED coupling constant at q 2 = m? z where only the quark and 
lepton contributions to the deviation from its q 2 = value, a(0) = 1/137.036, are included. Its 
magnitude is commonly referred to in the literatures |^9[ [30| It is denoted as a(m 2 z )f in 



cfs. 1 23, g5|, and is related to the coupling a(m z ) that contains the IF-boson contribution by 



l/a(m 2 z ) = l/a(m z z ) + 0.15 
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m t (GeV) |3 


173.8 ± 5.2 
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0.119 ± 0.002 
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128.90 ± 0.09 
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Table 1: Electroweak measurements at LEP, SLC and Tevatron. The average W- 
boson mass is found in ref. 32]. The reference SM predictions and the corresponding 
'pulP factors are given for m t = 175 GeV, niH SM = 100 GeV, ot s {m z ) = 0.118 and 
l/a(m z ) = 128.90. Correlation matrix elements of the Z line-shape parameters and 
those for the heavy-quark parameters are found in ref. 111]]. The data R e and vlpg are 
obtained by assuming the e-[i-T universality^ and are not used in our x 2 analysis. 



predictions of the MSSM are then expressed as the sum of the SM predictions and 
the difference between the predictions of the SM and those of the MSSM. This 
separation is useful because the SM predictions include parts of the known two- 
and three-loop corrections whereas the MSSM contributions are evaluated strictly 
in the one-loop order in the present analysis. 



2.1 Observables at Z-pole experiments 

The amplitude for the decay process Z — > f a f a is written as 

T{Z^fJ- a )=M[e z -J fct , (2.2) 

where e z is the polarization vector of the Z-boson and Jj a = fofl^fa = fl^Paf is 
the fermion current with definite chirality, a = L or R. The pseudo-observables of 
the Z-pole experiments are expressed in terms of the scalar amplitudes with 
the following normalization ]TT| : 

" kV2G F ml 0.74070' 1 ' ) 



A convenient parametrization of the effective couplings in generic electroweak the- 
ories has been given in refs. [24, 25]: 



s/2 = 


0.50214 + 0.453 Af z 
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(2.4a) 
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(2.4h) 



Here the first numerical terms in the r.h.s. of the above equations are the SM 
predictions at the reference point, (m t (GeV), m# SM (GeV), l/a(m z ), a s {m z )) = 
(175, 100, 128.90, 0.118), Ag z and As 2 are the universal gauge-boson-propagator 
corrections p3| , and Ag[, (a = L,R) are the process specific corrections. In the 
SM the following universality relations hold very accurately: 

(2.5a) 
(2.5b) 
(2.5c) 
(2.5d) 
(2.5e) 
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SM- 
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Only the term (AgfJsM has non-trivial m t and mn SM dependence. In the MSSM, 
we find that all the {Ag^) terms are non-vanishing, and the flavor-universality 
holds only among the first two generations. We study Ag%~ , Ag~£, Ag T R , Ag b L and 
Ag b R separately in the following sections. 

The universal part of the corrections Ag z and As 2 are defined as the shift in 
the effective couplings g z (m 2 z ) and s 2 (m|) []23| from their SM reference values at 
(m t (GeV),m HsM (GeV),l/a(m|)) = (175, 100, 128.90): 

f z (m 2 z ) = 0.55635 + Ag 2 z , (2.6a) 
s 2 (m|) = 0.23035 + As 2 . (2.6b) 

We find it convenient to express the above shifts in the two effective couplings in 
terms of the two parameters ASz and ATz, 

Ag 2 z = 0.00412AT Z , (2.7a) 
As 2 = 0.00360AS Z - 0.00241AT Z . (2.7b) 



The parameters Sz and T z are related to the S and T parameters [p 



S z = S + R - 0.064x a , (2.8a) 

T z = T + lA9R-—£, (2.8b) 

a 



in the notation of refs. p3| , p4| p5| . A compact summary of the definitions of 
the effective charges of ref. [23| and the oblique parameters S, T, U, R are given in 
Appendix B. The shift R 

An Arc 1 1 

-^—--^— = — R= — (1.1879 + AR), (2.9) 

accounts for the difference between the T parameter which measures the neutral 
current strength at the zero-momentum transfer ]19[ [33| and the Tz parameter 
which measures the quantity on the Z-pole. The x Q -dependent term in Sz reflects 
the fact that only a combination of S and l/a(m 2 z ) is constrained by the Z-pole 
asymmetry experiments p3| , p5| , p4]| . The factor 5 G denotes the vertex and box 
corrections to the muon decay constant 

GF = *V2m 2 w ' (2 ' 10) 

which takes 5 G = 0.0055 in the SM. In the MSSM, vertex and box corrections 
affect its magnitude and we define the shift AS G as 

~6 G = 0.00550 + A~§q. (2.11) 
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Although the supersymmetric contribution to A5 G is the correction specific to the 
muon decay rather than a part of the universal gauge-boson-propagator correc- 
tions, it appears in all the predictions of the electroweak observables because we 
use the muon decay constant Gp as one of the basic inputs of our calculations. 
We therefore include the effect A5 G as a part of our "universal" parameter, Tz- 

from the 



By inserting 



into 



and by using the shifts AS and AT p 



reference SM values of the S and T parameters, we can express the parameters 
A Sz and AT Z as 



AS Z = S z - 0.955 = AS + AR- 0.064a;, 



ATz = T, 



2.65 



AT +1A9AR- 

a 



(2.12a) 
(2.12b) 



The SM contribution to the shift AS*, AT and Ai? are parametrized as [23, 24 



(AS) S M 

(AT)sm 



-0.007x t + 0.091x h - 0.0104, 

(0.130 - 0.003x h )x t + 0.003x t 2 - 0.079^ - 0.028^ 



-0.0026a4 



(A^sm = -0.124 log 



26 



m HsM (GeV) 



log 



/ 26 

V100J 



(2.13a) 

(2.13b) 
(2.13c) 



and we set (A<5 G )g M = in our calculation. The explicit form of the MSSM 
contributions to the parameters AS", AT, Ai? are given in Appendix [B|. The MSSM 
contributions to the muon-decay parameter A5 G is given in Appendix [D], and those 
to the effective Z-boson decay amplitudes, Agf, are given in Appendix [Cj. 

The pseudo-observables of the Z-pole experiments are then obtained by using 
the above 12 effective couplings gl; 8 couplings of eq. ( |2.4j ) that are distinct in the 
SM, and the 4 additional coupling, g^J , g T L) g T R and g R , which can be distinct in 
the MSSM. The partial widths of the Z-boson are 



r = G F m% f 
1 3V2tt I 



2 C 



fv 



+ 



/ / 
9l~9r 




(2.14) 



where the factors Cfy and CfA account for the final state mass and QCD correc- 
tions for quarks. Their numerical values are listed in Table [| The c^-dependence 
in C q v,C q A is parametrized in terms of the parameter x s ( |2.1d| ). The last term 
proportional to a(m%)/ir in eq. ( 2.14|) accounts for the final state QED correction. 
The total decay width Tz and the hadronic decay width Th are given in terms of 
the partial width Tf: 



T z = 3T U + T e + T fl + T T + T h , 
T h = T u + T c + T d + T s + T b . 



(2.15a) 
(2.15b) 
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Table 2: Numerical values of factors Cfv, CfA for quarks and leptons used in eq. ( 2 . 



The ratios Ri, R c , R b and the hadronic peak cross section cr° are given by: 



Ri 



Rr 



Rb 



i2tt r P r, 



m\ T z 



2 ' 



(2.16) 



where I = e, \i or r. 

The left-right asymmetry parameter A* is also given in terms of the effective 
couplings gf, as 



A f 



{gif 



{g f R ? 



(g[) 2 + (gk) 



f\2 ' 



(2.17) 



The forward-backward (FB) asymmetry A^ B and the left-right (LR) asymmetry 
A° LR are then expressed as follows: 



^FB 



A 



LR 



A e . 



(2.18a) 
(2.18b) 



The effective parameter sin 2 measured from the jet-charge FB asymmetry is 
defined as 



sin 



2 o i t 



i g e R 



2g R -gi 



(2.19) 



All the Z-boson parameters in Table [T] are now calculable for arbitrary values 
of m t ,mH SM ,C(('m 2 z ) and a s (m z ), or x t ,x h ,x a and x s , respectively, in the SM by 
using the parametrizations given in this section. The predictions of the MSSM are 
calculated by using the formulae in Appendices [FJ, and |D|, by using the mixing 
matrix and the coupling notation [F28J of Appendix A. 
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2.2 The VF-boson mass 

The theoretical prediction of m w can be parametrized as |23], 

m w {GeV) = 80.402 + Am^, 
Am w {GeV) = -0.288AS + 0.418AT + 0.337A[/ 



+ 0.012a; - 0.126- 



A8, 



G 



a 



(2.20a) 



(2.20b) 



Here in addition to the S and T parameters, the U parameter |TP| is needed to 
calculate the effective charge p^(0) [EBJ that determines the muon decay constant 
Gf\ see eq. ( j2.10j ). The SM contribution to the shift Am w from the reference 
prediction can be parametrized as 125 



(Am iy ) SM = 0.064x t - 0.060z h - 0.009x| + 0.001x t (x t 
+0.001x| + 0.012x a , 



(2.21) 



which is obtained by inserting the SM contributions to AS, AT and AU [24, 25 



to eq. (|2.20b|) . The MSSM gives additional contributions to S, T, U and to A5 G : 

AS, 



Am 



w 



(Am 



vkJsm 



0.288S ncw + 0.418T ncw + 0.337f/ n 



0.126- 



G 



a: 



,(2.22) 



where 

S n 



AS - (AS)sm, T new = AT — (AT)sm, U ncvr = AU - (AC/) SM .(2.23) 



The MSSM contribution to S ncw , T new , C/ new and A5 G are given in appendices |B| 
and 0. 

2.3 Constraints on m t ,a(m 2 z ) and a s (m z ), and the SM fit 

All the electroweak observables are now calculated in the SM as functions of the 
four parameters, m t1 mH SM1 ot{m 2 z ) and a s (m z ), or x t ,x h ,x a and x s , respectively 
via ( |2.1D . The Higgs boson has not been found yet and the lower mass bound 



m HsM (GeV)^90 



(x h > -0.11) 



(2.24) 



is obtained from the LEP2 experiment . 

In our analysis we use the following constraints on the parameters m t 
a s {m z ) ]2{| and a(m z ) 



m t (GeV) = 173.8 ± 5.2 (x t = -0.12 ± 0.52), 
l/a(m|) = 128.90 ±0.09 (x a = ± 1), 
a s (m z ) = 0.119 ± 0.002 (x s = 0.33 ± 0.67), 



(2.25a) 
(2.25b) 
(2.25c) 



10 



as shown in the bottom of Table 0. As an alternative to the model-independent 
estimate ( |2.25b| ) f30|| , we also examine the case with the estimate ||31|| : 



l/a(m|) = 128.94 ± 0.04 (x a = 0.44 ± 0.44), (2.26) 

which is obtained by using the perturbative QCD constraints down to the r-lepton 
mass scale. 

We find that the reference point of our analysis (x t = x h = x a = x s = 0) is 
not far from the global minimum of the SM fit. With the external constraints of 
eq. (TO , we find ^/(do.f.) = 18.2/(22 - 4) at 



mi (GeV) = 172.3 ±5.0, 

a s (m z ) = 0.119 ±0.002, 
l/a(m|) = 128.90 ±0.09, 



m^ SM (GeV) = 1171 



98 
64- 



(2.27a) 
(2.27b) 
(2.27c) 
(2.27d) 



By replacing the estimate (P^5B|) |5p] by ( ggBJ ) |TJ, we find Xmin/( d -°- f •) 
18.3/(22-4) at 



m t (GeV) = 172.6 ±4.9, 

a s (m z ) = 0.119 ±0.002, 
l/a(m%) = 128.94 ±0.04, 



™# SM (GeV) 



14411. 



(2.28a) 
(2.28b) 
(2.28c) 
(2.28d) 



In summary, the SM gives a good fit to all the electroweak data if the Higgs boson 
mass is relatively light, with the mass below a few hundred GeV, as suggested by 
the ranges (|2727d|) or ( ggfcj ). 



3 The minimal supersymmetric standard model 
(MSSM) 

In subsection 3.1 we briefly summarize the MSSM and our coupling conventions. 
In subsection 3.2, we summarize the constraints on supersymmetric particle masses 
from direct search experiments. 



3.1 The MSSM Lagrangian 

We study consequences of the MSSM under the following constraints. 

(i) It has the minimal particle content, with the gauginos of the SU(3)c, SU(2) £ 
and U(l)y groups, three generations of squarks and sleptons, and two pairs 
of Higgs doublets and their superpartners. 
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(ii) It has the minimal superpotential, which is sufficient to give masses to all 
the quarks and leptons. In particular, we do not consider i?-parity non- 
conserving interactions. 

(hi) The three gaugino masses, M 3 ,M 2 ,Mi are taken to be independent, while 
most of our numerical results are obtained under the 'unification' condition, 
M 3 /M 2 /Mi = a 3 (m z )/a2{m z )/ai(m z ), where are the MS couplings. 

(iv) The scalar masses are introduced in such a way that squark and slepton 
interactions with neutralinos and gluinos are flavor conserving in the basis 
where quarks and leptons have definite mass. 

(v) The scalar masses of the first two generations are taken to be equal. We 
therefore have 5 scalar masses for the first two generations, and another 5 
for the third generation squarks and sleptons. 

(vi) We neglect mixings between the left- and right- chirality sfermions in the 
first two generations. Accordingly, we retain the soft SUSY breaking A 
parameters only in the third generations, A t , Af, and A T . 

Under the above constraints, the MSSM interactions can be parameterized in terms 
of 19 parameters: the ratio of the two v.e.v.'s tan/3, the Higgs- mixing mass /x, the 
pseudo-scalar Higgs-boson mass m^, the three gaugino masses M 3 , M 2 , Mi, the 
three A parameters A t ,Ab,A T , the five sfermion masses for the first two gener- 
ations, mg, rug, mg, m^, mg, and the five sfermion for the third generations, 
m Qz ' m u 3 > m r> 3 ' 171 l 3 ' 171 e 3 ■ Summing up, the 19 parameters of our MSSM La- 
grangian are 



All our analytic expressions are valid when the above 19 parameters of the MSSM 
are independently varied. Because the parameter space of the MSSM is too large 
even with the above restrictions, we present our numerical results often by vary- 
ing only the few most relevant parameters while keeping the rest of the parame- 
ters fixed at some appropriate values. Systematic investigation of the parameter 
space of the two representative models, the supergravity mediated and the gauge- 
interaction mediated supersymmetry breaking models, will be reported elsewhere 



tan/?, /i, rriA, 
M 3 ,M 2 ,M l , 
At, A},, A T , 



Q3 ,m U 3 > m D 3 ' m L 3 ' m E 3 - 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 
(3.1e) 



37|. 



12 



The MSSM Lagrangian with the above constraints are adopted for the new 
Feynman amplitude generator MadGraph2 |36 |, and its explicit form can be found 



in ref. P8[ . A compact summary of our notation is given in Appendix A. All the 
physical masses and couplings of the supersymmetric particles are calculated in 
the tree-level by using the MSSM Lagrangian with the above restrictions. 

We denote the two chargino mass eigenstates as, \i with m~- < m~- , the four 
neutralinos as, Xi > with m~o < • • • < m~o, and the gluino as g. The 7 sfermion 
mass eigenstates for each generation are denoted by, ul, d^, ur, cIr, V h l L , l R for the 
first two generations, and ti, t 2 ,bi,b 2 , v T , fx, r 2 for the third generation. The four 
Higgs boson mass eigenstates, the light and heavy CP-even neutral Higgs bosons 
h and H {m h < m H ), the CP-odd neutral Higgs boson A, and the charged Higgs 
boson H ± are obtained by using the improved effective potential of ref. [12] that 
assumes CP invariance[|. 

A few comments are in order. Although we do not consider models with non- 
minimal interactions, such as the .R-parity violating models or models with gauge 
interaction mediated supersymmetry breaking where the very light gravitino be- 
comes the lightest supersymmetric particle (LSP), all our results should be valid in 
those models with additional interactions, as long as their strengths are negligibly 
small as compared to the gauge interactions. Our analytic expressions allow for 
arbitrary CP violating phases in the mass parameters, [i and Aj, and hence some 
of the couplings have associated complex phases. Our numerical results are, how- 
ever, given in the CP conserving limit of the MSSM interactions. This is mainly 



because we adopt the effective Higgs potential of ref. [|12|] that assumes CP invari- 
ance in the scalar sector. Effects of CP violating interactions in the MSSM will be 
reported elsewhere. 

3.2 Bounds on the masses of SUSY particles from the di- 
rect search 

The predictions of the MSSM for all the electroweak observables reduce to those 



of the SM with the restricted Higgs boson mass range, <, 135 GeV O, Wi 



IB] , TE ], in the limit where all the supersymmetry breaking mass parameters are 
large. In the following analysis, we explicitly demonstrate this decoupling behavior 
quantitatively for all the electroweak observables. 

The differences between the predictions of the SM and those of the MSSM are 
hence largest when the supersymmetric particle masses are near their present lower 



t The radiative corrections beyond the 1-loop level on the lightest Higgs boson mass have 
been discussed in refs. 15, [l6|]. Those effects shift the theoretical prediction on rrih in the 
leading order a few GeV. Since the result of our analysis, however, is not affected quantitatively 
by such a small shift of mh, we take into account the leading order correction on mt in our study 
for brevity. 
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bounds from the direct search experiments. We calculate the particle masses from 
the 19 parameters of the MSSM, eq. ( |3.1| ), and confront them with those lower 
mass bounds. 

Current limits on the masses of scalar leptons are obtained from the LEP2 
experiments as J37 

m~ R £ 88 GeV, (3.2a) 
m~ R £ 80 GeV, (3.2b) 
m~ R £ 69 GeV, (3.2c) 

where Tr has been assumed to be the mass eigenstate. The corresponding limits for 
the chirality-left sleptons are weaker and depend on the mass of the charginos. The 
lower mass bounds for scalar quarks and gluinos are obtained from the Tevatron 
search experiments as [[38[ 

m~ £ 212 GeV, (3.3a) 
m~ Z 173 GeV, (3.3b) 

when the squark masses are common for the 5 light flavors and both chiralities, 
and when either the squarks or the gluino are much heavier. The bounds on 
squark masses depend, however, on details of their mass spectrum and on their 
decay patterns. The (almost) model-independent lower mass bounds are found 
from the LEP2 experiments [ j37fl : 

m 7i Z 88 GeV, (3.4a) 
m~ bi ^ 76 GeV. (3.4b) 

For charginos and neutralinos the following bounds are found from the LEP2 
experiments [pTjj : 

m~o Z 33 GeV, (3.5a) 

m~- £ 90 GeV. (3.5b) 

Those on the Higgs particles are |p5 



m H SM > 95 GeV, (3.6a) 

m h ;> 84 GeV, (3.6b) 

m A £ 85 GeV, (3.6c) 

m H - ^ 69 GeV. (3.6d) 
In addition, for the neutral particles, we assume 

m~ > 45 GeV, (3.7a) 

m~o + m~o > m z , (3.7b) 
so that they do not contribute to the total width of the Z-boson flTT|. 
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4 Oblique corrections in the MSSM 



As explained in the introduction, the MSSM contributions to the vertex or box 
corrections vanish whenever either the sypersymmetric scalars (squarks and slep- 
tons) or the supersymmetric fermions (charginos, neutralinos, and the gluino) are 
heavy enough. In those cases, the MSSM particles can still affect the electroweak 
observables via their contributions to the gauge-boson propagators, which are of- 
ten called the oblique corrections |0| |2(J, [22f . Because the oblique corrections 
affect all the electroweak observables in a flavor-independent manner (universal- 
ity), and because their effects are found to be most significant under the present 
constraints on the new particle masses as summarized in the previous section, we 
study them in this section in great detail. 

The formalism presented in section 2 tells us that the precision electroweak 
experiments at the Z-boson pole constrain just two oblique parameters, Sz and 
Tz, whereas the P^-boson mass m w can be taken as the third oblique parameter. 
We favor m w over the U parameter as our third oblique parameter, because we 
can avoid correlations among the three oblique parameters this way and because 
we could not gain insight by adopting the U parameter in our MSSM analysis. 

In the following subsections, we examine the oblique corrections from each 
sector of the MSSM, since all the oblique corrections are expressed as linear sum 
of individual contributions [^6|: squarks and sleptons (sec. 4.1), the MSSM Higgs 
bosons (sec. 4.2), and charginos and neutralinos (sec. 4.3). The effects of combining 
all the contributions are discussed in subsection 4.4. As remarked above, the 
dominance of the oblique contributions to the electroweak observables is justified 
only when either the sfermions or the supersymmetric (-ino) fermions are heavy. 
All the oblique contributions are of course relevant always as a part of the full 
MSSM contributions. 

In presenting our results, we find it most convenient to parameterize first the 
experimental constraints on the oblique parameters, and then confront the MSSM 
contributions to those parameters against the experimental constraints. We further 
observe that the uncertainty in the SM contributions to the oblique parameters 
which arise from the uncertainty in rriH SM Q2.24| ), m t (|2.25a ) and a(m|), (|2.25b| ) or 
( |2.26| ), is significant as compared to the magnitudes of the MSSM contributions. 
The uncertainty in a s {m z ) ( |2.25c|) has been found to affect our results little. We 
therefore obtain the constraints on the oblique parameters in such a way that we 
can examine the sum of the SM and the MSSM contributions for arbitrary values 
of m t ,m HsM and a{m 2 z ) within their present constraints. This can be achieved by 
observing that the m t and m# SM dependences of the SM predictions appear only in 
the oblique parameters and in the Zb L b L vertex correction, Ag b L , and by observing 
that their a(m 2 z ) dependences appear only through the combination Sz ( 2.12a ) 
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Figure 1: Constraints on ASz — 33.7 Ag b L and ATz — Q0.3Ag b L from the electroweak 
precision measurements. The symbol (x) denotes the best fit from the electroweak 
data. The 39% (Ax 2 = 1) and 90% (A% 2 = 4.61) contours are shown. The SM 
predictions are given for m t = 165 ~ 180 GeV and m# SM = 90 ~ 150 GeV. 

and Am w ( gggg ). 

By adopting the 5 parameters, ASz, ATz, Ag b L , Am w and a s {m z ) as the free 
adjustable parameters, we find the following fit for all the electroweak observables 
of Table 1 , under the constraint (|2.25c|) : 



p = 0.89, (4.1a) 



AS Z - 33.7 Ag\ = -0.070 ± 0.113 1 
AT z -60.3Ag b L = -0.183 ± 0.137 J 

Am w ( GeV) = 0.008 ± 0.046, (4.1b) 



2 iC , /A^ + 0.00086\ 2 
Xmm = 15.4 + , (4.1c) 



V 0.00076 

where d.o.f.= 20 — 5 = 15. 

The results of the fit are shown in Figs. [I] and 0, where the SM predictions for 
some representative values of m t and rriH SM are given for l/a(m%) = 128.90 (x a = 
0) p0[ . The predictions for different estimate of l/a(m|) can easily be obtained 
by using the ^^-dependences of ASz (|2.12a|) and Am w ( |2.20b| ) . Fig. ^] shows 



the constraint ( fTTa| ) in the plane of AS Z - 33.7 Ag b L and AT Z - m.3Ag b L . It 
is clearly seen from this figure that relatively low value of m t as compared to 
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Figure 2: The total x 2 versus m w . The 1-er experimental bound on m w is 
shown. The SM predictions are given for m t = 165 ~ 180 GeV and rriH SM = 100 ~ 
150 GeV. The reference SM point at m t = 175 GeV, rriH SM = 100 GeV, l/a(m|) = 
128.90 is shown by the arrow. 

the present experimental mean value (|2.25a|) and relatively light Higgs boson are 
favored within the SM fit. In Fig. [| the SM prediction for m w (at x a = 0) and the 
resulting total \ 2 are shown. The m w data slightly favors rnn SM = 100 GeV over 
m H SM = 150 GeV for m t <, 170 GeV, so does the total x 2 . When m t ^ 175 GeV, 
the trend is reversed and m# SM = 150 GeV is favored against m^ SM = 100 GeV in 
the total x 2 -Q 

The SM fits reported in section 2, eqs. fl2.27p and ( |2.28| ), are obtained from the 
above fit by combining it with the constraint on m t ( |2.25a| ), and that on l/a(m|), 
( ^.25b| ) or ( |2.26| ), respectively. In both Figs. ^ and^|, the theoretical predictions are 
given for the reference value of l/a(m|) = 128.90, or x a = 0, in eq. (|2.1d ). Effects 
of changing l/a(m 2 z ) can easily be studied by shifting the theoretical predictions 

* The Xtot value in the figures ||, |i| |^ for the SM reference point at m t = 175 GeV and 
m H SM = 100 GeV is about 0.25 larger than the quoted value of 19.8 in Table |l|, because of the 
contribution from the a s (m z ) constraint (|2.25c). 
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shown for x a = by 



AS z (x a = 0) -> AS z (x Q = 0) - 0.064x a , (4.2a) 
Am w/ (x Q = 0) Am w (x a = 0) + 0.012a; Q , (4.2b) 



as required by eqs. fl2.12"aD and ( |2.20b| ). Change in l/a{m z ) of the full uncertainty 



of the conservative estimate Q2.25bp makes x a = ±1 in eq. ( |2 . 1 c| ) . We can clearly 
see from Fig. [I] that the resulting horizontal shift of the SM prediction due to 
eq. (|4.2aj) can affect significantly the preferred range of the Higgs boson mass in 



the SM. Inspection of Fig. § with the shift due to eq. (|4.2b|) tells us that the effect 
is not significant for m w with its present measurement error. We can verify the 
effect by comparing the SM fits (|2.27| ) and ([2.281) , that were obtained by using the 



two estimates, (|2.25b ) and ( p. 26 ), respectively. The mean value of the estimate 
( |2.26 ) is about x a = 0.4, and hence the SM prediction in Fig. [I] moves horizontally 
in the negative direction by about 0.03. This may affect the Higgs boson mass 
by about 30%. Although naive, this simple estimate reproduces qualitatively the 
difference in the most favored value of m HsM in the two fits, that are found to be 
about 22% between eqs.(^27|) and (|2~2^ ). 



In the following subsections, we show the MSSM contributions to the oblique 
parameters by superposing them on Figs. [I] and 0, by choosing the reference SM 
point at m t = 175 GeV and m# SM = 100 GeV (x t — x h — 0). The MSSM 
predictions for other choices of m t and mH SM are then obtained simply by shifting 
the SM reference point. 

4.1 Squarks and sleptons 

In Fig. |3] we show the sfermion contributions to ASz and ATz that are superposed 
to the SM contribution of Fig. [l|. The origin of the plot (ASz — ATz — Ag L = 0), 
marked by the big open circle, gives the SM prediction at m t = 175 GeV and 
m HsM = 100 GeV for l/a(m 2 z ) = 128.90 (x a = 0). The contributions of the 
squarks and sleptons are shown separately. The net result should be obtained by 
adding the SM contribution, the squark contribution and the slepton contribution 
vectorially in the two-dimensional plot. 

It is clear from the figure that for m t ~ 175 GeV and mn SM ^ 130 GeV, both 
the squark and slepton contributions make the fit worse than the SM. The situation 
does not improve by changing our estimate for l/a(m' z ) in the range (|2.25b|) . The 



squark contribution makes the fit worse, because it always makes Tz larger than 
the SM prediction, which is already larger than the preferred value of the data. 
The slepton contribution makes the fit worse, because it gives negative Sz with 
either T z = (tan/? = 2) or with slightly positive T z (tan/3 = 50). The better 
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Figure 3: The squark and slepton contributions to (ASz, ATz) for tan/3 = 2 
and 50. The SUSY breaking scalar masses for the left-handed and^right-handed 
squarks are assumed to be same, denoted by m^. The ti-tj? and &l-&r mixings 
are controlled by A cS = A* ff = A b cS . In the slepton sector, although sizable f L -f R 
mixing may be induced for large tan f3, the mixing is neglected because it does not 
affect both ASz and ATz significantly. The reference SM point (ASz = ATz = 
Ag\ = 0) is marked by the open circle at the origin. 

fit to the data requires a contribution with both ASz < and ATz < 0, which 
cannot be achieved by the contributions from squarks and sleptons. 

Let us examine the squark and slepton contributions to the ASz and ATz 
parameters in more detail. We find that they do not contribute significantly to R, 
and hence we can understand the qualitative behavior by studying their contri- 
bution to S and T. To fix our notation, we first give the squared mass matrix of 
the squarks and sleptons. Since the structure of the mass matrix of sfermions are 
similar among the different flavors, we give the mass-squared matrix of the stop 
as an example. In the (tL,tn) basis it is given by 

M| = 1 mX Ijf ' ' (4 ' 38) 



m~ L = m-^ + m| cos2/3(I 3u - s 2 Q u ) + m 2 t , (4.3b) 
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m~ 



m~ + m| cos 2(3s 2 Q u + m 2 



u 3 



t ! 



(4.3c) 



where s 2 = sin 2 #vk(^z)ms- The diagonal elements of the mass-squared matrices 



for the other sfermions are obtained by replacing 



for u, 



{m§,m 3 



for d, 



{m l3 ,m^ 



for r, 



(m z , m~ 



by 



for 6, 



for e, and by choosing 



appropriate values for the third component of the weak isospin I^f and the electric 
charge Qf. The off-diagonal elements in ([4.3a]) , m t v4* ff , should be replaced by 

m f Af s 



A f 



Af - fi* cot p {f = t) 
Af — /j,* tan (3 (f = b,r) 



(4.4) 



for b and f, while it is set to zero for the other sfermions. The sfermion mass- 
squared matrix can be diagonalized by using the unitary matrix U* for f — t,b, f: 



(U f YMfU f 



diag(m? 



m~ < my , 
y h h' 



(4.5) 



In all our numerical examples, we put real values for Af and fi. Since we neglect 
the left-right mixing in the first two generations, their mass eigenvalues are given 
by the diagonal elements. 

It is useful to examine analytic expressions of the S and T parameters under 
some approximations. If the mixing between the left- and right-handed states is 
negligible, the sfermion contribution to the S parameter, AS", can be given by the 
following simple form: 



AS 



1 



nm 2 z 



T, C f^f Y fL F 5(m 2 z : mj L , 
f 



m 



(4.6) 



where Cf is 3 for the squarks and 1 for the sleptons. The symbol / runs over 
the flavor space. In this limit, only the SU(2)l doublets contributes to the S 
parameter. The explicit form of the function F 5 (m 2 z : mj, m~) in the first line has 
been given in ref. p3]| . Since AS* is proportional to the hypercharge Yf, the relative 
sign of AS is opposite between squarks (Y q = 1/6) and sleptons (Yl = —1/2) in 
the first two generations. Moreover, since AS is proportional to I%f , it receives 
contribution only when there is a mass splitting among SU(2)l multiplet members. 
It takes a particularly simple form in the large sfermion mass limit. For example, 
the ul-cLl and Vi-Il contributions can be expressed as 



UL 



m~ 



12n m~ + m~ 

UL d L 

i mA — mA r 
12tt mA" + mi- 



O 



1 + 



1 1 

mA — m~ 

U L d L 

mi- + m~ 

U L d L 
1 1 

mA — mA 
v i II 



mA 



2 

mA , 
II 



(4.7a) 
(4.7b) 
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respectively. The magnitude of AS is determined by the mass difference between 
the up- and down-type states in the SU(2)i doublet. Except for the ti-bi case, 
the mass difference of the up- and down-type components of the SU(2) doublets 
is determined by the D-terms in the mass-squared matrix by neglecting the tiny 
contribution from the fermion masses 

m~ — m~ = m~ — m~ = (1 — s 2 )m% cos 23. (4.8) 

Since the r.h.s. of eq. ( |4.8|) is negative for tan/3 > 1 and its magnitude grows 
as tan 3 increases, we can understand the qualitative behavior of the squark and 
slepton contributions in Fig. [3|. 

The analytic form of the T parameter in the zero left-right mixing limit of 
squarks and sleptons are also simple. It also receives contribution from particles 
which carry the SU(2)x quantum number. For example, the t^-b^ contribution is 
given by 



V2 47r 2 a ' 



2 2 2 

mA mA mA 
2 i ™2 \ , ft b L , b L 



-(mi • ni~ ) + ,: L "\ In 



2 *i b L mi — mi mi 

th b L t L 

F q v t L b L 



G ^(rr^-miy 



24\/27r 2 a" / mi +rr£ 

tL b L 



1+0 



2 2 

mA — mA 
t L b L 

mi + mi 



(4.9) 



tL b L 

The second line in eq. ( |4.9| ) is an expression in the large tL, Sl mass limit. Contri- 
butions from the other squarks and sleptons are obtained by replacing the masses 
and Cf factors appropriately. The squarks in the first two generations and sleptons 
contributions can be found by replacing the mass parameters and Cf appropri- 
ately. Since the function in the parenthesis of eq. (|4.9|) is positive for any values 
of {m~ L , mj L ) or (m~ , mj^ ) , AT always receives positive contributions from the 
squarks and sleptons. 

The special attention should be payed for the stop-sbottom contribution. There 
is a large mass difference between tL and b\, which is dominated by m 2 — m 2 . It 
is rather easy to understand that the effect of the left-right mixing of the stop 
and sbottom qualitatively. Since the off-diagonal elements of the mass-squared 
matrices of the stop and sbottom are multiplied by the top and bottom mass 
squared, respectively, the effect of the left-right mixing in the stop sector is much 
larger than that in the sbottom sector. Let us focus on the left-right mixing in the 
stop mass matrix for brevity. Suppose that only t\ is light and the other squarks 
are heavy and degenerate. Suppose also that there is no mixing between 6^ and 
With these assumptions, the squark contribution to AT in the third generation 
may be parametrized by and m^ (~ mj 2 w m^ ) and the unitary matrix U t . 
Then, the stop-sbottom contribution to AT is give by 



AT w ° F 1 C 
y/2 4ir 2 a ' 



\U T n \ 2 F 5 (0 : m Ti ,m~ ) - I^H^I^O : m^m~ ) 
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Figure 4: The squark and slepton contributions to m w for tan (3 = 2 and 50. The 
SUSY breaking scalar masses for the left- and right-handed squarks/sleptons are 
taken to be common, which are denoted by msusY- Two cases of the left-right 
mixing, A e s = GeV (left) and A e s = 300 GeV (right) are examined. In each 
figure, two cases of the SM reference point, (m t (GeV), m# SM (GeV)) = (175,100) 
and (170, 100), are marked by circles. 
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(4.10) 



Let us recall that the left-right mixing is induced by A l eE in eq. (|4.4j) . The factor 
| Ulil 4 decreases (< 1) as A l eS increases and it suppresses AT. The behavior of 
the squark contribution in Fig. |3| may be interpreted in this way. On the other 
hand, if only is kept small in eq. ( 4.10|) while keeping \U\-\_\ finite, the T 
parameter grows with m~ . This reflects the growth of the Al s parameter as 

(rr^/mJKlVa 



It/^l 2 ), that breaks the custodial SU(2) symmetry 



Let us give the analytic expression of AR in the large sfermion mass limit. 
Because AR is a linear combination of the Z-boson propagator corrections between 
two different momentum transfer scales, there are contributions not only from the 
left-handed sfermions but also from the right-handed sfermions. By taking mj^ 
and mj R as the left- and right-handed sfermion masses, respectively, AR in the 
heavy mass limit is given by 



AR 



30tt 



-C t 



(% - s 2 QfY-f + (s 2 Qf) 



2 m l 



ra- 



ni' 



fa 



(4.11) 
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Although the negative sign of AR leads to negative contribution to ATz in ( |2.12b|) , 
its magnitude is so tiny that the net contribution of sfermions to ATz is found to 
be always positive. 

In Fig. [I], we superimposed the sfermion contribution to m w and to the total x 2 ■ 
A e ff = and A e g = 300 GeV cases are shown side by side. Because the total x 2 is 
not a linear sum of the SM contribution and the new physics contribution, we show 
two representative cases of the reference SM point, one for (m t , m HsM ) = (175, 100) 
and the other for (m t ,m HsM ) = (170, 100) in GeV units. For each representative 
(m t ,mH SM ) case, we show the MSSM predictions of m w and the total x 2 as a 
function of the common sfermion mass msusY- Here, for simplicity, we set all 
the 10 sfermion mass parameters in eqs. (|3.1d|) and (|3.1e|) to have a common value 



t^susy, and the 3 A( s parameters in eq. (|4.4j ) to have a common value A e g. We find 
that the sfermion contributions always make m w larger than the SM prediction. 
Because the SM prediction for m w is smaller than the experimental mean value at 
m t 170 GeV, the sfermion contribution can improve the fit for m w . As may be 
seen from the four examples in Fig. |], however, the sfermion contribution always 
make the overall fit worse than the SM. By comparing the A e $ = case (left) and 
the A e g = 300 GeV case (right), we find that the unfavorable sfermion contribution 
can be made small by introducing A e s for the same value of mgusY- The trend 
can be understood qualitatively by using the analytic expressions (|4.7|), ( |4.9| ) and 
( |4. 1Q ) for the sfermion contributions to the 5* and T parameters. The sfermion 
contribution to the U parameter is found numerically small. 

As is clear from the figures, the present m w data gives little contribution to 
the total x 2 °f the fit. Although the positive contribution to Am w can improve 
the fit to the m w data if m t ~ 170 GeV, the overall x 2 always gets worse because 
of the Z-parameter constraints as summarized in Fig. [3[ 

4.2 MSSM Higgs bosons 

The MSSM has an extended Higgs boson sector with three neutral Higgs bosons, 
h, H, and A, and one charged Higgs boson H . All their masses and couplings are 
calculated in terms of the MSSM Lagrangian parameters, by using the improved 
one-loop potential |Tj|. Their contribution to the oblique parameters are evalu- 



ated carefully, because parts of the SM corrections contain the m# SM dependence 
from the two-loop corrections. They appear in the T parameter and in the Ag L 
parameter, although the m HsM dependence of the latter turned out to be negligibly 
small. In order to avoid discontinuity in the theoretical predictions because of the 
lack of the corresponding two-loop results in the MSSM, we make the following 
arrangement in our actual calculation. For each MSSM parameters, we first cal- 
culate the lightest CP-even neutral Higgs boson mass, m^. We then evaluate all 
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the SM corrections for m# SM = m^. The contributions of the MSSM Higgs sector 
to the radiative corrections are then evaluated as the sum of the SM contribu- 
tions, that partially contain the two-loop correction, and the difference between 
the MSSM and the SM Higgs sector contributions that is calculated strictly in the 
one-loop order by using the formulae in the Appendices. In this way, we can test 
numerically the decoupling of the MSSM Higgs sector in the large limit. 

In Fig. 5, we show the contributions of the MSSM Higgs sector to the ASz 
and ATz parameters when itla = 50 GeV, 100 GeV and 300 GeV. As always, the 
results for tan (3 = 2 are shown by solid blobs, and those for tan (3 = 50 are shown 
by open blobs. We find that the MSSM prediction is remarkably near to the SM 
prediction at mj SM = rah already at rriA = 300 GeV. The lightest Higgs boson 
mass in ref. |12] can be approximately expressed as 
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Figure 5: Contributions to (ASz, ATz) from the Higgs bosons for tan (3 = 2 and 50. 
The 1-loop improved Higgs potential, in which only the top and bottom Yukawa 
couplings are retained, are used to find the mass eigenvalues. The inputs are the 
pseudo-scalar Higgs-boson mass tua and tan/5. Three values of tua, ttia — 50, 100 
and 300 GeV are examined. The scale parameter msusY which appears in the 
1-loop potential is set at 1 TeV. 
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18.8 
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0.023 


-0.015 


-0.014 


18.7 




1000 


129 


1000 


1003 


0.025 


-0.017 


-0.015 


18.7 



Table 3: Oblique parameters, ASz, AT Z and Am w in the MSSM Higgs sector for 
tan (3 = 2 and 50. The mass spectrum of the Higgs bosons are calculated by using 
the 1-loop improved scalar potential which is approximated by retaining only the 
top- and bottom-quark Yukawa couplings |12|] . The scale parameter m SU sY which 
appears in the 1-loop potential is set at 1 TeV. 



+20x 



SUSY 



1-0.08o; SUS y 1 + 1.5 1- 



tan (3 



(4.12) 



when rriA = 300 GeV. The parameter xsusy is defined by xsusy = hi(msusY/l TeV). 
Eq. (|4TT2| ) is valid for 2 < tan/3 < 50 and 1 TeV < m SU SY < 2 TeV in which the 
error is smaller than 0.6 GeV. 

We find m h = 102 GeV and 129 GeV for tan/3 = 2 and 50, respectively, at 
rriA = 300 GeV for m SU sY = 1 TeV, and the predictions are remarkably near to 
the SM predictions at m# SM = m^. The rriA dependence of the MSSM predictions 
are shown in Table 3, together with the masses of all the MSSM Higgs bosons for 
tan (3 = 2 and 50 and at msusY — 1 TeV. From the Table, we can see that the 
predictions of the MSSM reduces essentially to those of the SM at rnn SM = rn h 
when m A £ 200 GeV. 



4.3 Supersymmetric fermions: charginos and neutralinos 

In this subsection, we study the contributions of the charginos and neutralinos. 
The chargino and neutralino mass matrices depend on four parameters, the two 
gaugino masses Mi and M 2 , the supersymmetric Higgs- mixing mass \i and tan/3. 
Because the gaugino masses and the Higgs-mixing mass are invariant under the 
electroweak gauge symmetry, they contribute neither to the S nor the T parameter. 

On the other hand, because they are fermions, their virtual creation can affect 
the running of the gauge-boson propagator, the R parameter ( |2.9| ), strongly when 
the pair-creation threshold is near the Z-boson pole |23], [39]. 
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Figure 6: The chargino and neutralino contributions to (ASz, ATz) for tan/3 = 2 
and 50. The supersymmetric Higgs-mixing mass \i is fixed at fi — 300 GeV. The 
GUT relation for the gaugino masses, M 2 /&2 = M\jo.\ is assumed, and cases 
for M.2 = 140 ~ 90 GeV are shown by using the marked SM reference point at 
(m t , Tn,H SM ) = (175,100) in GeV unit as the origin. Those points with the cross 
( x ) symbols give the lighter chargino mass m~- below 90 GeV. 

Although the present lower mass bound on the charginos is as large as the 
Z-boson mass rather than its half, we find that the effect can still be significant. 
We show in Fig. |6| the contributions of the charginos and neutralinos to the ASz 
and ATz parameters. In the figure, the supersymmetric Higgs mass is fixed at 
\x = 300 GeV, and the SU(2) L gaugino mass M 2 has been varied. The U(l)y 
gaugino mass is scaled by the relation Mi/M? = cti/&2 for definiteness. The case 
for tan/? = 2 are shown by solid blobs, and those for tan/5 = 50 are shown by 
open blobs. As the M 2 decreases, the lightest chargino mass decreases, and those 
points with the cross symbol give the chargino mass below the direct search bound 

(HI). 

Most importantly, we find that the effect of relatively light charginos and neu- 
tralinos gives both ASz and ATz negative, which improves the fit over the SM. 
This is essentially because of the R parameter that reside both in the ASz and 
the ATz parameters ( |2.12j ). For instance, if the wino-like chargino mass is near to 
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Figure 7: The chargino and neutralino contributions to (rn w ,Xtot) f° r tan/3 = 2 
and 50. The supersymmetric Higgs-mixing mass /i is fixed by // = 300 GeV. The 
GUT relation for the gaugino masses, M 2 /cb = M\ja.\ is assumed. The 6 blobs 
are for M 2 (GeV) = 140,130,120,110,100,90 at tan/5 = 2 (solid blobs) and at 
tan/3 = 50 (open blobs). Those points with the cross (x) symbols indicate that 
the lightest chargino mass m~- is smaller that 90 GeV. 



half the Z-boson mass, its contribution to the parameter R behaves as 

'1 16 
j§ ~ 37 



Aflw-c 4 ^-^), (4.13) 



where c 2 is given by c 2 = 1 — s 2 . [3 = JiM$/m% — 1 is the analytic continuation 
of the chargino velocity below the threshold. When 4M|/m| 3> 1, the effect is 
suppressed as 

4c 4 mi 



AR ps — — — ^ 
15tt M 2 



(4.14) 



'2 

The chargino contribution to AR is negative just like the sfermion contribution 
to AR in eq. Q4.ll ). The coefficient of the wino contribution in eq. ([4.14 ) is 



found to be about 90 times larger than that of the Ir contribution. This large 
negative contribution to AR makes both ASz and ATz significantly negative 
when a relatively light chargino of mass about 100 GeV exists. 
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Table 4: Chargino and neutralino contributions to the oblique parameters, 
AR, ASz, ATz, Am w and the total x 2 °f t ne fit at tan (3 = 2. 

More detailed numerical results are summarized in Table 4 for tan/? = 2 
and in Table 5 for tan/5 = 50. From both tables, we can see that ATz is 
always negative for these inputs and its magnitude is dominated by AR, see 
eq. ( f2.12b|) . On the other hand, ASz and Am w can be positive or negative. As 
compared to the total x 2 value of 20.1 which is obtained at the SM reference value 
(m t (GeV), mH SM (GeV), a(m|) = (175,100,129.90) under the a s {m z ) constraint 
( j2.25c| ), see eq. ( |4.1| ), we find an improvement of Xtot m the tables by between 1.0 
and 1.9, by including the contribution of a chargino if its mass is around 100 GeV, 
just above the present mass bound. 

4.4 Summary of oblique corrections 

We show in Fig. |8| the contributions of all the sectors of the MSSM together. The 
total oblique corrections are obtained by summing the SM contribution for one's 
favorite m t at rriH SM = 100 GeV, the MSSM Higgs sector contribution, the squark 
contribution, the slepton contribution, and the chargino-neutralino contribution 
vectorially in the plane. 
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Table 5: Chargino and neutralino contributions to the oblique parameters, 
AR, AS Z , AT Z , Am w and the total \ 2 of the fit at tan /3 = 50. 

It is surprising that each sector of the MSSM gives distinctive contributions 
in the ASz and ATz plane. The squark sector contributes essentially to positive 
Tz direction, which is strongly disfavored by the electroweak data. The slepton 
sector contributes negatively to Sz, but Tz remains constant or slightly positive 
for large tan f3. This is also found to be disfavored by the data. Therefore, if only 
the squarks and sleptons are light, their contributions would make the fit to the 
electroweak data significantly worse than the SM. 

On the other hand, the MSSM Higgs sector is found to give the contribution 
very similar to that of the SM when the lightest CP-even Higgs-boson mass m h is 
taken to be the SM Higgs boson mass. Therefore, if rrih is significantly larger than 
the reference value of 100 GeV, either with large tan /3, or by having very heavy 
squark masses that affect the effective scalar potential, the fit improves slightly. 

Finally, the contribution of charginos and neutralinos make both ASz and ATz 
negative, and hence the fit improves more significantly if the lightest chargino mass 
is near the present direct search limit of around 100 GeV. 

As emphasized in section 2, the oblique corrections dominate the whole ra- 
diative effects if either squarks and sleptons, or the charginos and neutralinos are 
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Figure 8: Summary of the supersymmetric contributions to (ASz, AT^). 

very heavy. Otherwise, in addition to the oblique corrections summarized in this 
Figure, there will be additional vertex and box corrections. 



5 Non-oblique corrections 

We found in the previous section that the existence of relatively light charginos 
(of mass ~ 100 GeV) improves the SM fit to the electroweak data, whereas the 
existence of relatively light squarks and sleptons (of mass ~ 100 GeV) makes 
the fit significantly worse. On the other hand, if both supersymmetric fermions 
(charginos and neutralinos) and supersymmetric scalars (squarks and sleptons) are 
light, then in addition to the universal gauge-boson-propagator corrections (the 
oblique corrections), we should expect process specific vertex and box corrections 
to be non- negligible. We would like to address this problem in this section. 

There are just two types of process specific corrections which are relevant for 
the precision electroweak experiments on the Z-boson properties and the P^-boson 
mass. One is the vertex corrections to various Z — > ff amplitudes, the parameters 
Ag^, in eqs.( |2.4|) . Even if we assume the universality of the first two generation 
squarks and sleptons, there are 14 independent Zff vertex corrections in the 
MSSM. The other is the vertex and box corrections to the muon-decay amplitude, 
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the parameter A5 G , that affect the ATz parameter ( f2.12b|) in the Z-boson exper- 
iments and the Am w parameter ( |2.20b| ) for the H^-boson mass. Although this is 
just a single correction factor, it affects all the electroweak predictions because we 
use the observed magnitude of the muon decay constant as one of the basic 
inputs of our calculation. 

We present our studies on these non-oblique corrections step by step in the 
increasing order of the number of observables that will be affected. In section p.l\ 
we study the MSSM Higgs boson contribution to the Zbb, Ztt and Zv T v T vertices. 
In section |5.2j , we study the order a s vertex correction when gluinos and squarks 
are both light. This will affect the Zqq vertices only, and we find it useful to study 
its quantitative significance independently of the other electroweak corrections. 
In section we study the electroweak corrections to the Zqq vertices when 
both squarks and charginos/neutralinos are light. Effects on the Zbb vertices are 
carefully studied. In section |5.4j , we study the Zll vertices when both sleptons and 
charginos/neutralinos are light. We notice, however, if both the supersymmetric 
fermions (charginos and neutralinos) and the sleptons are light enough to affect 
the Zll vertices significantly, they should also affect the muon-decay amplitude, 
A5 G , and hence all the remaining electroweak observables. All the numerical 
results in the following subsections are found at (m t (GeV), a s (m z ), a{m z )) = 
(175,0.118,128.90). 



5.1 Z/ /-vert ex corrections by the MSSM Higgs bosons 

When tan (3 is very large, the Yukawa couplings of the 6-quark and the r-lepton get 
large, and we expect significant vertex corrections for both Zbb and Ztt vertices. 
The most significant contributions are found to arise from the vertices with the 
Higgs bosons. In Fig. |, we show Rb, T T and the total x 2 as functions of m^, 
the pseudo-scalar Higgs-boson mass in our Higgs potential. The tan/5 = 2 case 
is given by the solid lines, whereas the tan/? = 50 case is given by the dashed 
lines. For simplicity, we set all squarks and slepton masses as well as the lightest 
chargino mass to be 1 TeV so that their contribution can be neglected. 

We find that for sufficiently small uia at large tan/3, there is a significant 
positive contribution to both and T T . These effects have been noted |^0[ when 
the experimental data on Rb larger than the SM prediction was reported ||41|| . We 



find that the fit to the latest Rb data is slightly better with larger Rb predicted in 
the MSSM with tan j3 = 50 and ~ 60 GeV, but as we can see from the total 
X 2 behavior, the scenario is not favored over the SM. This is because the MSSM 
predicts light Higgs particles whose contributions to the oblique parameters worsen 
the fit, as shown in Table 3 in Section || Also the pseudo-scalar Higgs-boson mass 
ttla significantly below 100 GeV is not acceptable because it gives the lightest 
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Figure 9: The Higgs boson contributions to Rb (top), T T (middle) and the total \ 2 
(bottom) as functions of itla for tan/? = 2 (solid line) and 50 (dashed line). The 
squarks, sleptons and lighter chargino masses are fixed by 1 TeV. The other parameters, 
m t ,a s (m z ) and a{m 2 z ) are fixed by 175 GeV, 0.118 and 128.90, respectively. The 1-a 



allowed range of Rb and T T (83.94 ± 0.22 MeV [11]) are shown in the top and middle 
graphs, respectively. Two thick lines in the bottom graph denote the SM fit for rriH SM = 
106 GeV(solid) and rriH SM = 129 GeV(dashed) which are the lightest Higgs boson mass 
in the MSSM at via = 1 TeV for tan/3 = 2 and 50, respectively. 



m 



CP-even Higgs boson mass (m/J below the direct search bound fl3.6|). 

In Fig. |] we show by horizontal lines the Xtot °f the SM prediction at 
175 GeV and l/a{m%) = 128.90 for m HsM = 106 GeV and 129 GeV. These Higgs 
boson masses correspond to the masses of the lightest CP-even Higgs boson (h) 
in the MSSM at tan/3 = 2 and 50, respectively. Table 3 shows that the oblique 
corrections from the MSSM Higgs sector reduces to the SM predictions at the 
corresponding Higgs boson mass, m# SM = m h , already at uia = 300 GeV. In con- 
trast, Fig. p| shows that the vertex corrections can affect the total x 2 significantly 
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Figure 10: The SUSY-QCD correction to Re as a function of the gluino mass 
M3 = M~. All the squark masses are assumed to be given by the universal scalar 

maSS, TTlsquark- 

even at rriA = 300 GeV, especially for tan j3 = 50. This is mainly because of the 
enhanced bbH and ttH couplings at large tan/?, where H is the heavier of the 
CP-even Higgs bosons with mass m H « m A . 

5.2 SUSY-QCD corrections 

Supersymmetric QCD corrections to the Z-boson hadronic width, I\, or its ratio 
Rg to the leptonic width Ti is one of the first vertex corrections to the Z-boson 
decay amplitudes calculated ||. They contribute at the order a s level, and hence 
could be significant if both squarks and gluinos are light. We would like to first 
examine their quantitative significance in view of the present lower mass bounds 
on squarks and gluinos fl3.3|). 

The one-loop Zff amplitudes of Appendix C, which is presented by using a 
generic notation that apply for all the vertices and for all the MSSM corrections, 
take a particularly simple form when only the order a s corrections are considered. 
For instance, if the mixing between the left and right chirality squarks is negligible, 
which we assume to be the case for the first two family squarks and sleptons, we 
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Inputs (GeV) 


0(a s ) vertex 


vertex + 


obi. 




m u L rn~ u M~ 


R e /R e (SM) 


XsQCD 


Rt/R e (SM) 


Xtot 


I 


100 100 100 


1.00152 


24.1 


1.00241 


129.1 


II 


180 180 180 


1.00048 


20.5 


1.00095 


52.9 


III 


1000 180 180 


1.00003 


19.9 


1.00005 


20.5 


IV 


1000 180 421 


1.00001 


19.8 


1.00004 


20.4 



Table 6: The SUSY-QCD corrections to R e /R e (SM) and x 2 for tan/3 = 2. In 
all cases, m~ fl = mj is assumed. The universality of the SUSY breaking scalar 
masses among the different generations, = ttIq s , = m ij 3 anc ^ m d = 171 d 3 * s 
also assumed. The left-right mixings in the stop and sbottom sector are taken to 
be zero by using an appropriate values of A t , A b and \i. The effects of the SUSY- 
QCD vertex correction for four cases are summarized in the second column while 
the SUSY-QCD vertex correction and the oblique correction by squarks are given 
in the last column. Only in IV, the GUT relation M~ = M3 = (a^/ ' dt-^M^ has been 
respected. The value of M~ in IV leads to M 2 = 120 GeV, and m~- = 94 GeV for 
/i = 300 GeV. 



find 



z 

2C 24 (m! : m~ ,M~,m~ , 



(5.1) 



where the two- and three-point functions, B , B\ and C24 are given in ref. [53|. We 
reproduced the results of refs. M. 

In Fig. |IU| we show the SUSY-QCD correction to the ratio Ri as functions of 
the gluino mass M 3 = M~ when all the squark masses are taken to be the same for 
brevity. Significantly larger magnitude of Re as compared to the SM prediction at 
the same a s is found only when the common squark mass is as light as 50 GeV and 
the gluino is lighter than 100 GeV, consistent with the early observation || . Once 
we take into account the present lower mass bound of squarks and gluinos fl3.3|) , 
the quantitative significance of the SUSY-QCD correction diminishes. 

We show in Table || the magnitude of the SUSY-QCD correction to Re for 
several squark and gluino mass inputs. We should compare the percentage correc- 
tions with the present experimental accuracy of Re, which is 0.125% from Table |l| 
Only when all the squarks and gluino masses are around 100 GeV, we can ex- 
pect non-negligible effect from this sector. When all the masses are moved up to 
180 GeV, consistent with the present direct search limits ( |3.3|) , the effect on Re 
is about 0.05% and is already insignificant. Our studies on the oblique correction 
suggest that the left-handed squarks should be rather heavy not to worsen the 
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electroweak fit. If only the right-handed squarks are kept light at 180 GeV, to- 
gether with the gluino of the same mass, the effect essentially goes away, as shown 
for the case III in Table |6|. In the last line of the Table, we study the case IV 
where M~ = M3 = 421 GeV, which is the mass obtained by using the 'unification' 
relation M 3 /M 2 = in one of our selected MSSM points that accommodate 

a chargino of mass around 100 GeV. The effect diminishes to 10~ 4 . 

In the following analysis, we assume that the left-handed squarks are relatively 
heavy (ra~ J> 1 TeV) in order not to spoil the good fit to the electroweak data. 
We also adopt the 'unification' relation for the three gaugino masses M3/&3 = 
M2/&2 = Mi/di in all the subsequent numerical examples. Under such conditions, 
we find that the SUSY-QCD corrections cannot have significant effect on our 
electroweak studies. 



5.3 Electroweak Zqq vertex corrections with squarks 

In this subsection we explore the possibility that the electroweak SUSY corrections 
to the Zqq vertices improve the overall fit to the data. As in the case of the SUSY- 
QCD corrections with a gluino-squark loop, we find that m~ L « mj should be kept 
high ( ^ 1 TeV) in order for their positive T contribution not to spoil the good 
overall fit of the SM. 

The only case that SUSY electroweak correction to the Zqq vertices are found 
to give nontrivial effects is when loops of light charginos and a light £1 affect the 
Zbb vertices. We show in Fig. [IT] the effect of supersymmetric vertex corrections 
to the ratio R^. In order to avoid large unfavorable oblique corrections to the T 
parameter, we set all the q~L masses large by setting = = 1 TeV. We 
take the singlet squark masses at = = = mg 3 = 200 GeV near the 
Tevatron search limit ( |3.3| ). The lighter stop, £1, can still be light enough to affect 
the Zb L b L vertex by having large £l-£r mixing. 

In the figure we show the present l-a experimental bound on R^, and the 
MSSM predictions as functions of the £1 mass. The lower part of the figure shows 
the total x 2 of the fit that includes contributions from all the electroweak data in 
Table [I]. Three cases with a light chargino (m~- <^ 100 GeV) have been examined: 

(i) Xi is almost wino [p = 300 GeV, M 2 = 120 GeV, m~- = 94 GeV) 

(ii) Xi is almost higgsino (ji = 100 GeV, M 2 = 800 GeV, m~- = 93 GeV) 

(iii) Xi is almost higgsino (pi = 90 GeV, M 2 = 200 GeV, m~- = 54 GeV ) 

The case (iii) violates the LEP chargino mass bound ( |3.5| ). As is clear from these 
figures, the MSSM effects on the Zbb vertices are not significant under the present 
constraints on the £1 and Xi masses. The scale of the total x 2 -> which is significantly 
below that of the SM, is set by the light chargino-neutralino contribution to the 
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Figure 11: The stop contribution to Rb (top) and the total \ 2 (bottom) as 
functions of for = 1 TeV and = = 200 GeV where the universality 
in the generation space of the squark sector is assumed. Three different inputs for 
the chargino mass are studied: (i) (ji,M 2 ) = (300,120), (ii) (//, M 2 ) = (100,800), 
(iii) (//, M 2 ) = (90, 200) in the GeV unit. All the other SUSY particle masses are 
fixed by 1 TeV. The SM prediction, Rb = 0.21566 and x 2 — 19.8, are also shown. 

oblique parameters as we explained in the last section, and the slight decrease that 
we observe for the cases (i) and (ii) for lighter vnq is not significant. Only when 
both the t\ and Xi are as light as 50 GeV and when the Xi is almost higgsino and 
t\ is almost t R , we could obtain significant positive contribution to Rb as discussed 
in the literatures a couple of years ago [S, [10| . We reproduced all the numerical 



results presented in ref. [10 



With the above choice of squark parameters, = 1 TeV and mjj = = 
200 GeV, and with our neglect of all the left-right mixing in the first two genera- 
tions, the electroweak vertex corrections to the four light-quark flavors are found 
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Figure 12: The MSSM contribution to A5 G /a (a) and the total \ 2 (b) function 
of the left-handed slepton mass, mj for tan (3 = 2. In the figures, fi and M2 are fixed 
at 300 GeV and 120 GeV, respectively. The GUT relation for the gaugino masses are 
assumed. In (a), the vertex and box corrections are shown by dotted and dashed lines, 
respectively. The total contribution, 2(5 vcr tcx + 3box, is given by the solid line. In (b), 
the right-handed slepton mass is fixed by mj = 100 GeV. The masses of squarks and 
Higgs bosons besides the lightest Higgs are fixed by 1 TeV. The solid and dashed lines 
show the total \ 2 with and without AS G , respectively. 

to be negligibly small. 

5.4 Z//-vertex and ^-decay corrections with sleptons 

So far, we set all the slepton masses to infinity so that they do not participate in 
the MSSM corrections. When sleptons are relatively light, together with relatively 
light charginos and neutralinos as motivated by their favorable contributions to 
the oblique parameters ASz and ATz, we expect two distinctive effects in the 
electroweak predictions of the MSSM. First, they affect the muon-decay amplitude 
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Figure 13: The scalar- lepton contributions to Re (top), A e (middle) and Ap B (bottom) 
as functions of m~ for tan [3 = 2. In each figure, the impact of the right-handed slepton 

'■L 

mass is shown by solid (mj = 100 GeV) and dashed (m~ = 1 TeV) lines, respectively. 
The other parameters are the same with those in Fig. [l2|. 

and hence the parameter AS G /a, in eq. ( j2.10| ), that appear in the basic expression 
for the oblique parameter ATz ( |2.12b|) and the IF-boson mass (|2.20b|) . Second, 
they affect the Zll vertices, and hence the parameters, Ag l L , Ag l R , Ag u ^ in eqs.( [2.4|) . 

We first examine the MSSM contributions to the muon-decay amplitude care- 
fully, because if they can make A5 G /a positive ATz can also become negative: see 
eq. fl2.12b|) . A combination of this negative ATz from the muon-decay correction 
and the negative ASz from the oblique corrections (see Fig. 3) could improve the 
overall fit to the Z-pole data. If the muon-decay amplitude receives positive contri- 
bution from new physics, then it is equivalent to make the effective charged current 
strength larger, and the T parameter, which is a measure of the neutral-current 
to charged-current strengths, is effectively reduced. On the figures of the oblique 
parameters ASz vs. ATz, the positive A5 G /a would be regarded as the negative 
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contribution to the ATz parameter, which is favored by the electroweak data. On 
the other hand, as soon as the muon-decay amplitude is affected, we should expect 
non-oblique Zll vertex corrections to affect the MSSM predictions for the Z-boson 
parameters. We therefore present in Fig. 12 the MSSM contribution to A5 G /a 
and that to the total x 2 , as functions of the II mass by assuming = ^j 3 - 

By using the analytic expressions for the muon-decay amplitude in Appendix 
C, we find that box corrections are always positive, while the vertex corrections 
are negative. After summing up, we find that the net MSSM correction to A5 G 
can become positive only when mr < 130 GeV, when the lighter chargino mass 

<>L/ 

is around 90 GeV. Fig. [l| shows that the total % 2 improves slightly by including 
the A5 G /a correction to the muon decay, Xtot (withA5 G ) < Xtot (without A5 G ), 
when mj ^130 GeV and A5 G < 0. 

Fig. [13] shows the slepton contributions to R^, A e , A^ B when tan/5 = 2 and the 
lighter chargino mass is around 90 GeV (/i = 300 GeV, M<i = 120 GeV). In order 
to find the impact of the right-handed slepton on the observables, we show the 
cases of mr = 100 GeV and mj = 1 TeV by solid and dotted lines, respectively. 
We find that the non-oblique corrections from the slepton-ino loops are generally 
small and that they cannot remedy the unfavorable oblique effects of light left- 
handed sleptons. On the other hand, the right-handed slepton mass mj is not 
constrained significantly by the electroweak data. 



6 Summary and Discussions 

In this paper, we studied radiative corrections to the electroweak observables due 
to new particles in the MSSM systematically, confronting them against the lat- 
est data on the Z-boson parameters and the IV-boson mass. After introducing 
our analytic framework and giving parametrizations of the SM predictions, we 
first studied the effects of the oblique (gauge-boson-propagator) corrections from 
squarks and sleptons (sec. 4.1), the MSSM Higgs bosons (sec. 4.2), charginos and 
neutralinos (sec. 4.3), separately. The effects of non-oblique corrections are then 
studied systematically; the MSSM Higgs-boson contributions to the Zbb and Ztt 
vertices (sec. |5.1|) , the gluino-squark-loop corrections to the Zqq vertices (sec. 15.21) , 
the squark-chargino/neutralino-loop corrections to the Zqq vertices (sec. [5.3D , and 
the slepton-chargino/neutralino-loop corrections to the Zll vertices and the muon- 
decay amplitude (sec. |5.4| ). When studying contribution of each sector of the 
MSSM, we set all the irrelevant new particle masses sufficiently high so that their 
contributions can be neglected. 

Our observations can be summarized as follows. The agreement between the 
SM predictions, and the electroweak data of Table 1 is excellent: for instance, 
the four-parameter fit ( |2.27D to the 22 data points of Table |TJ (18 Z-parameters, 
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m w ,m t ,a s (m z ) and a{m z ) P0|) finds x^ in /(d.o.f.) = 18.2/(22 — 4), with the 



probability 44%. However, close examination of the fit shows that the data favor 
additional small negative contributions to the two oblique parameters, Sz and Tz- 
see Fig. |l[ Small positive contribution to m w can also improve the fit slightly if 



< 175 GeV: see Fig. f§ 



Our study of the oblique corrections in the MSSM can be summarized as follows 
by using, as a measure of the goodness of the MSSM fit as compared to the SM 
fit, the quantity 

A X 2 = (xL)mssm - (Xtot)sM, (6.1) 

where we compare the two fits by using the common parameter set (m t (GeV), 
l/a(m|), a s {m z ) ) = (175, 128.90, 0.118) at m HsM =m h = 100 GeV. 

• Small left-handed slepton mass {Tnj) makes ASz negative while keeping ATz 
essentially unchanged. It contribute positively to Am w . The fit worsens 
slightly (Ax 2 ^l) if m-^ 170 GeV for tan/3 = 2 or mj L ^250 GeV for 
tan/3 = 50, and significantly (Ax 2 ^4) if mj <>110 GeV for tan/5 = 2 of 
m T ^150 GeV for tan/? = 50. 



• Small left-handed squark mass parameters (mg, m^) make ATz and Am w 
positive while keeping ASz essentially unchanged. The effective t L -t R mixing 
parameter A\ s of the magnitude of the order of mg 3 tame the unfavorable 
positive ATz contribution, but not completely. The fit worsens significantly 
(Ax 2 ^ 4) if m^(= m^ 3 ) ^ 340 GeV for tan/3 = 2 or m^(= m^) ^ 300 GeV 
for tan (3 = 50, even for A cS = 300 GeV, which correspond to 300 GeV 
(tan/3 = 2) and m 7i ^ 260 GeV (tan/3 = 50). 

• The MSSM Higgs boson contribution can be approximated by the SM Higgs 
boson contribution at m HsM = m h if the pseudo-scalar Higgs-boson mass 
satisfies tua ^ 200 GeV. The effects on the electroweak observables are small, 
| Ax 2 1 ^ 1 for tan/3 = 2 or |Ax 2 | ^0.5 for tan/3 = 50, when the direct search 
bound of rrih > 75 GeV (|3.6c|) is taken into account. See Fig. 5 and Table 3. 

• In contrast, the light charginos and neutralinos are found to contribute neg- 
atively to both ASz and ATz, while they contribute negligibly to Am w . 
Their contribution can hence improve the SM fit (Ax 2 ~ —1): see Figs. ^| 
and |7|. 

We find that the best fit is obtained when the lightest chargino mass (m~- ) is near 
its experimental lower bound, m~- ~ 90 GeV, see Tables 4 and 5. Overall picture 
of the oblique corrections in the MSSM is summarized in Fig. RL 
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Non-oblique corrections in the MSSM are then studied systematically in sec- 
tion H as follows. The MSSM Higgs bosons can affect the Zbb, Ztt and Zv T v T ver- 
tices through their Yukawa couplings to quarks and leptons of the third generation 
(sec. |5.1| ). The SUSY-QCD corrections to the Zqq vertices are studied when both 
squarks and gluinos are light (sec. |5.2| ). The contributions of the light squarks and 
light charginos/neutralinos to the Zqq vertices are then studied (sec. pTB| ). When 
both sleptons and charginos/neutralinos are light, not only the Zll vertices but 
also the muon-decay amplitude is affected (sec. We find: 



The vertex corrections due to the MSSM Higgs bosons can be significant 
(A% 2 ^ 4) if the pseudo-scalar Higgs-boson mass uia is smaller than about 
100 GeV for tan/? ~ 1. For tan/3 ~ 50, the effects are significant when 
m A ^ 300 GeV. The overall fit to the electroweak observables can only be 
worse than the SM: see Fig. |9|. 

The order a s SUSY-QCD corrections to the Zqq vertices are found to be 
negligibly small (|Ax 2 |<^0.5) when squarks and gluino masses are bigger 
than about 200 GeV: see Fig. [10] and Table. |6|. 

The Zqq vertex corrections due to light squarks and light charginos/neutralinos 
are found to be insignificant if we take into account the direct search mass 
bounds on the lightest squarks ( |3.4j ) and the lightest chargino (|3.5|). We also 
found that good overall fit to the electroweak data with m~- ~ 100 GeV can 
be maintained with m-jr ~ 100 GeV if the left-handed squark mass is kept 
large, ^q(= m Q 3 ) ^ 1 TeV: see Fig. [LI. 



• When the left-handed slepton mass is small (mj <^ 200 GeV) and charginos/neutralinos 
are light (m-- ~ 100 GeV), the Zll vertices as well as the muon-decay am- 
plitude are affected significantly. Although we find that their contribution 
to the muon-decay amplitude can improve the fit when mj ~ 100 GeV, its 
unfavorable oblique contributions make the overall fit worse than the SM 
(Ax 2 £ 1) for mj L ^ 150 GeV. 

Summing up, we find that no contribution of the MSSM particles can improve 
the fit to the electroweak data by taking into account the non-oblique correc- 
tions to the Zff vertices and the muon-decay amplitude. Best fit to the data 
is found when the supersymmetric fermions (charginos and neutralinos) are light, 
m~- ~ 100 GeV, and the five scalar mass parameters mA,rn^,m^,m^ and m^ 3 
are all large; ^ 300 GeV, m~ = ^ 500 GeV, = m~ g ^ 1 TeV. The 
singlet scalar-mass parameters, m^^m^m^m^m^ and do not affect the 
electroweak fit significantly, and hence light SU(2)i-singlet scalar bosons are al- 
lowed. 
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Before closing, we give a few examples of the MSSM predictions for all the 
19 electroweak observables of Table [I] when we have light charginos/neutralinos 
(m~- ~ 100 GeV) and relatively light singlet squarks and sleptons. We examine 
the four cases I to IV of Table [7], all of which have a light chargino of about 
100 GeV and the lightest squark at around 200 GeV, and the lightest slepton 
at around 150 GeV. For brevity, we assumed that the SUSY breaking scalar 
masses, 777^,771^,771^,771^,771^ are universal among different generations, and the 
effective left-right mixing mass parameters of t,b and f, A 1 ^, A b eS , A^ s , are also 
common. The doublet squarks and sleptons are kept heavy at around 1 TeV and 
the pseudo-scalar Higgs-boson mass rriA is set at around 500 GeV in all four cases. 
The scenarios I and II are for tan (3 = 2, III and IV are for tan f3 = 50. fi = 300 GeV 
for I and III while 11 = —300 GeV for II and IV. Masses of most relevant physical 
states are also shown in Table [7]. Because M 2 < in all four cases, the lightest 
charginos/neutralinos are dominantly gauginos. The lightest squarks and sleptons 
of the third generation squarks and sleptons all have predominantly right-handed 
component. The gluino mass (M~ = M 3 ) and the U(l)y gaugino mass Mi are 
fixed by the 'unification' condition M3/M2/M1 = a.3/0.2/0.1. 

The predictions are compared with the data in Table § for the scenarios I 
and II (tan/5 = 2) and in Table |^ for the scenarios III and IV (tan/3 = 50). In 
each Table we give the reference SM predictions at m t = 175 GeV and mn SM = 
100 GeV for comparison. The first column in each scenario shows the contributions 
of charginos/neutralinos and the MSSM Higgs bosons only. The results of this 
column represent the predictions of the MSSM when all new particles except the 
charginos/neutralinos and the Higgs bosons are heavy. In the second column, 
we show the predictions when only the oblique corrections from the squarks and 
sleptons are added, and in the last column we show the complete MSSM predictions 
including all the vertex and box corrections. We show in the second column 
the predictions of the oblique corrections only in order to show the quantitative 
importance of non-oblique corrections. 

From Table || we find that the MSSM predictions are almost completely deter- 
mined by the oblique contributions of charginos and neutralinos in the scenarios 
I and II with tan/? = 2, where the three scalar-mass parameters, ttla, ttl^ and 
are all kept large. The relatively light squarks and sleptons of a few hun- 
dred GeV barely contribute to oblique or non-oblique corrections. In contrast, the 
tan (3 = 50 cases of III and IV given in Table |9| show that the non-oblique correc- 
tions can have significant effects even though the physical mass spectrum in all 
four cases are similar in Table [7|. This is because the enhanced Yukawa coupling of 
the 6-quark allows the relatively light bi of predominantly right-handed component 
to contribute to the Zbb vertices. This affects r& (and hence Rb in Table |9|) and Th 
(and hence Re). Similar contributions to the Zt l tl vertex are found not to affect 
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Input 


s 












tan j3 


n 


M 2 


A t 


m A 


m Q 


m u 






m E 


T 
1 


2 


300 


125 


-657 


486 


956 


196 


181 


894 


129 


TT 
11 


2 


-300 


90 


568 


434 


1273 


183 


176 


639 


159 


TTT 
111 


50 


300 


110 


812 


481 


1091 


221 


204 


925 


165 


IV 


50 


-300 


105 


-657 


486 


956 


196 


181 


894 


129 




Outputs 




m h 




m~o 


M~ 


m~ 


m~ 


m 7! 


mr 

6i 


m u R 


mj 


I 


104 


98.2 


53.4 


439 


134 


893 


214 


182 


194 


182 


II 


109 


101 


48.3 


316 


162 


637 


231 


177 


180 


177 


III 


130 


101 


53.6 


386 


169 


922 


247 


196 


219 


206 


IV 


128 


98.4 


51.7 


369 


133 


892 


230 


168 


193 


183 



Table 7: Inputs and outputs parameters of the MSSM scenarios I to IV with light 
charginos/neutralinos and light squarks and sleptons. 



the fit significantly. 

Finally, we would like to perform the 'global' fit to all the electroweak data 
of Table [l] in the MSSM, by taking into account the present data on m t (|2.25a| ). 
l/a(m|) (|2.25b|) and a s {m z ) ( |2.25c| ). From the previous comprehensive studies, 
it is clear that the improvement of the total x 2 over the SM can be realized only 
when a light (~ 100 GeV) chargino exists and when the left-handed squark and 
slepton masses (m^ and mj\ and the pseudo-scalar Higgs-boson mass (m^) are 
all sufficiently large. In the following study, we assume that the universality of the 
scalar mass parameters among the different generations, mg = ^q 3 , m^} = m u 3 > 
nip = itl^ 3 , = m Z 3 an< ^ m E = 171 Ez~ Also, we fix all the scalar mass parameters 
at 1 TeV {wIq = m-^ = rn^j = rn^ = = itia = 1 TeV) and A{ s = for brevity. 
Then, there are three parameters left in the MSSM, tan /3, M 2 and /i. We find the 
total x 2 as a function of the lighter chargino mass m~- for tan/3 = 2 (Fig. [14]) 
and tan/3 = 50 (Fig. |i~5|) . Since it is worth studying how the mixing between 
the gaugino and the higgsino affects the fit, we show the following three cases 
separately: M 2 //i = 0.1 (solid lines), 1 (dotted lines) and 10 (dashed lines). The 
small number of M 2 //i implies that the lighter chargino is dominantly the wino 
while the large number of M 2 //i implies that it is dominantly the higgsino. The 
behavior of Xtot m ^ ne small m~- region ( <> 60 GeV) can be understood from the 
chargino contribution to the R parameter, eq. ( ^4.13| ), which behaves as 1/(3 when 
m~- is close to a half of m 7 . 

In our assumption on the mass spectrum of the MSSM, the lightest Higgs 
boson mass is predicted as nih = 106 GeV for tan/? = 2 and rrih = 129 GeV 
for tan/3 = 50. The decoupling in the large SUSY mass limit is examined by 
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Xtot 



19.8 



18.1 



18.2 



18.1 



18.2 



18.3 



Table 8: The pull factors of the predictions of the MSSM scenarios I and II. Those 
of the reference SM predictions (see Table. [|ll]]) are given for comparison. 
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SM* 



III 



inos/Higgs 



obi. 



Pull 



all 



IV 



inos/Higgs 



obi. 



.4 



Re 

FB 



A e 

Rb 
R c 
A o,b 

FB 
aO,c 
FB 



sin^f 



A 

An 



m w 

Xtot 



-1.4 
0.3 
0.7 
0.3 



f Re 


0.7 




1.3 


{ Rr 


-0.7 


f A% 
{ A ^ 


-0.5 
-0.1 
1.1 



-1.2 
-0.1 

1.2 
0.3 
-2.4 
-0.8 

0.7 

1.0 
-1.9 
-0.5 

0.2 
19.8 



-0.8 
0.3 
0.8 
0.6 

0.8 
1.3 
-0.6 

-0.4 
0.1 
1.2 

-0.9 
0.2 

1.2 
0.3 
-1.9 
-0.6 

0.6 

1.6 
-1.9 
-0.5 

0.3 
17.8 



-0.8 
0.3 
0.8 
0.6 

0.8 
1.3 
-0.6 

-0.4 
0.1 
1.2 

-0.9 
0.2 

1.2 
0.3 
-2.0 
-0.6 

0.6 

1.6 
-1.9 
-0.5 

0.3 
17.8 



-0.8 
0.2 
1.1 
0.6 

0.9 
1.5 
-0.5 

-0.4 
0.1 
1.2 

-0.8 
0.2 

1.4 
0.3 
-2.1 
-0.6 

0.6 

1.6 
-2.0 
-0.5 

0.3 
19.4 



-0.8 
0.3 
0.8 
0.6 

0.8 
1.3 
-0.6 

-0.4 
0.1 
1.2 

-0.9 
0.2 

1.2 
0.3 
-2.0 
-0.6 

0.6 

1.5 
-1.9 
-0.5 

0.3 
17.7 



-0.8 
0.3 
0.8 
0.6 

0.8 
1.3 
-0.6 

-0.4 
0.1 
1.2 

-0.9 
0.1 

1.2 
0.3 
-2.0 
-0.6 

0.6 

1.5 
-1.9 
-0.5 

0.2 
17.8 



Table 9: The pull factors of the predictions of the MSSM scenarios III and IV. 
Those of the reference SM predictions (see Table. ||1 1|| ) are given for comparison. 
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Figure 14: 



m~- for tan /3 
line. 



The total x 2 in the MSSM as a function of the lighter chargino mass 
2. The SM best fit (x 2 = 18.2) is shown by the dashed horizontal 

106 GeV 



The dot-dashed horizontal line shows the SM fit using m#, 



SM 



which is the lightest Higgs boson mass predicted in the MSSM. Three different 
M 2 -/i ratio (10, 1, 0.1) are studied. The bound on m~- from the LEP2 experiment 
is shown by the dashed vertical line. 

comparing the MSSM fit with the SM fit at ttlh sm = rrih rather than the SM best 
fit at m HsM = 117 GeV. We show by horizontal lines the corresponding SM results 
in Figs. Ill and |I5|. The decoupling seems to hold at m~- = 1 TeV for tan/? = 2. 

I 1 I 1 1 Xi 

On the other hand, the MSSM prediction at tan/3 = 50 differs slightly with the 
SM prediction at itlh sm = rrih even when m~- = 1 TeV, as seen from Fig. 15. We 



looked for the origin of this small discrepancy and found that it comes from the non- 
oblique corrections as shown in Table [TOI. Among the 14 distinct Zff vertices of 



the MSSM that appear in eq. (2J3), we find that the Zb R b R vertex receives the non- 
negligible corrections at tan/5 ~ 50. Its origin is the enhanced Yukawa coupling 
among the 6-quark, the t-quark and the heavy Higgs bosons (A, H, H ) of 1 TeV 
mass. In particular, the charged Higgs boson contributes to the Zb R b R vertex 
has an additional enhancement factor due to the exchanged top-quark mass. The 
Yukawa couplings among the charged Higgs boson, the right-handed 6-quark and 
the left-handed t-quark is proportional to m b tan ft, which gives the largest coupling 
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Figure 15: The total x 2 i n the MSSM as a function of the lighter chargino mass 
m~- for tan (5 = 50. The SM best fit (x 2 = 18.2) is shown by the dashed horizontal 
line. The dot-dashed horizontal line shows the SM fit using m HsM = 129 GeV 
which is the lightest Higgs boson mass predicted in the MSSM. Three different 
M 2 -fi ratio (10, 1, 0.1) are studied. The bound on m~- from the LEP2 experiment 
is shown by the dashed vertical line. 

in the MSSM when tan f3 > ^Jm t /mb ~ 7. The combined effects make the effective 
Zb R b R coupling smaller than the SM by about 0.6% even when m H - ~ 1 TeV. 
This 0.6% decrease in g b R slightly worsen the fit to the Rb and data. 

Taking account of the direct search limit on the chargino mass from LEP2, we 
can conclude that m~- ~ 100 GeV gives the best fit point of the MSSM in both 
small and large tan (3. Among the three ratios of M 2 //x, the gaugino dominant case 
(M 2 /fJ> = 0.1) gives slightly better fit at m~- ~ 100 GeV rather than the others. 
The improvement of the MSSM fit over the SM is not very significant, however, 
since the SM fit to the data is already good. 
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tan (3 = 2 


tan p = 50 




obi. 


obi. + vtx. 


obi. 


obi. + vtx. 


V 2 . 

A mm 


18.5 


18.5 


18.3 


18.7 


Xsm at m H SM = m h 




18.5 




18.3 



Table 10: Xmm * n the MSSM for tan/5 = 2 and 50 which are estimated by using 
the oblique corrections only (obi.) and all the MSSM corrections (obi. + vtx.). 
All the scalar mass parameters are the same with those in Figs. and 13. We 
take M2I '/J* = 0.1 as an example. As a reference, x 2 hi the SM (x|m) are shown 
for 77iHg M = 106 GeV and mn SM = 129 GeV, which should each correspond to the 
decoupling limit of the tan (3 = 2 and tan (3 = 50 cases, respectively. 



7 Conclusion 

We have presented the results of our comprehensive study of the constraints on 
the MSSM parameters from the electroweak precision measurements. The gauge- 
boson-propagator (oblique) corrections are parametrized in terms of three param- 
eters, ASz,ATz and Am w . The Sz parameter is obtained from the standard S 
parameter by taking into account the running effect of the Z-boson propagator 
correction between q 2 = and q 2 = m 2 z , which we denoted by R, while the T z 
parameter is obtained from the standard T parameter by taking account of the 
R parameter and A5 G , which is the correction to the muon-decay amplitude. We 
found that the left-handed sfermion contributions always make the fit worse than 
the SM due to the positive contributions to the T z parameter (squarks) or the 
negative contributions to the Sz parameter (sleptons). The contributions from 
the right-handed sfermions to the oblique parameters are negligible. The contri- 
butions of the MSSM Higgs bosons to the oblique parameters behave like the SM 
Higgs boson contribution at wh sm = mh when the pseudo-scalar Higgs-boson mass 
rriA is large (m^ ;> 300 GeV). The most important finding of our study is that the 
light chargino contribution can make the fit better than that of the SM due to its 
effect on the R parameter. We therefore studied vertex and box corrections that 
contain light charginos/neutralinos and find that no combination of squarks and 
sleptons can improve the fit further. 

The 'global' fit of the MSSM has been performed by assuming the heavy mass 
limit of all the sfermions and the pseudo-scalar Higgs boson. We found that the 
best fit is obtained when the chargino mass is at the current lower bound from the 
LEP2 experiments. 
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Note added 



After submission of this paper, we learned that Djouadi et al. |44]] calculated the 
two-loop QCD corrections to the gauge boson two-point functions. They found 
that the QCD correction enhances the one-loop squark contributions to the p- 
parameter (the T-parameter in our paper) by 10%-30%. This effect may shift 
effectively the stop mass in our results slightly higher. We would like to thank 
S. Heinemeyer for calling our attention to Ref. ||44[| . 

We also learned that it is not appropriate to use the experimental data on the 
effective weak mixing angle, sin 2 Q 1 ^ (see Table 1) which is measured from the 
jet-charge asymmetry, in our analysis when the non-universal vertex corrections 
are significant. This is because the parameter is obtained from the asymmetry 
data by assuming that the Z-boson decay branching fractions to each quark flavor 
obey the SM prediction. We confirm that our results presented in this paper do 
not change significantly when we remove the jet-charge asymmetry data in the fit. 
We are grateful to K. Moenig for clarifying the problem for us. 
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Appendix 



A Mass eigenstates and couplings in the MSSM 

All the analytic expressions for the 1-loop contributions of the MSSM particles are 
expressed in terms of the masses of the mass-eigenstates and their couplings. We 
summarize the notation of ref. [28] here for completeness. 

The mixing matrix elements between the current-eigenstates and the mass- 
eigenstates appear in the expression for the gauge-boson-propagator corrections. 
There are four types of the mixing matrices in our restricted MSSM. For squarks 
and sleptons, we consider only the chirality mixing among the third generation 
squarks and sleptons: 

h = E^Oii/i , fl = i:(U%f*, (A.la) 

3=1 3=1 

.//,• •%./}• r R = T,(u%f;, (A.ib) 

3=1 3=1 

for / = t, b and f. For the charginos we have 

2 2 

WE = E( U l)ijXJl , W+ = E( U lTiMr , (A.2a) 

3=1 3=1 
2 2 

H~ dL = E( U Lhxj L , H+ R = Y,{U2)%xU , (A.2b) 

3=1 3=1 
2 2 

Wr = J2(Ur)i3XJr , Wt = E(^)iiXi , (A.2c) 

3=1 3=1 
2 2 

H~n = EiUZhxjR > H+ L = J2(Ug)* 2j xt L , (A.2d) 

3=1 3=1 

and for the neutralinos we have 

Bl = Y,{U»)v% L , B R = )r(US)v% R . (A-3a) 

(A.3b) 

H" dL = E(^Wz , Kr = J2(U^h3xU , (A.3c) 

3=1 3=1 

H° uL = hu») 4j x% , H° uR = J2(U») 4j x% . (A.3d) 





)ljXjL > 




= E(^)i^, 


3=1 


































































)ijXjL > 






j'=i 
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Finally, the Higgs bosons are expressed as follows: 

H+ = (H-)* = i X + sin/5 + H + cos/3 , (A.4a) 
HJ = (H+)* = i X ~ cos/3 + H~ sin/3 , (A.4b) 

H® = —j= (v sin /3 + h cos a + H° sin a — ix 3 sin /3 + iA cos (3) , (A. 4c) 
V2 

H° d = -^= (v cos (3 - /i sin a + if cos a + ix 3 cos /? + iA sin /?) . (A.4d) 
V2 

Here the phases a and /3 are fixed by requiring cos /3, sin /3, cos a > 0. 

For notational compactness, we adopt the following generic notation of the 
MSSM couplings in Appendices [C] and |D] for the vertex and box corrections. Only 
three types of the couplings appear among the 9 types of the renormalizable cou- 
plings listed in ref. [p8| . The FFV (fermion-fermion- vector) couplings, g„ lF2V \ are 
defined from the corresponding Lagrangian term as 

C = g^ v T\rP a F 2 V fl (A.5) 

where a = L,R. The SSV (scalar-scalar-vector) couplings, g SlS2V , are defined by 

C = ig s ^ v S{d»S 2 V^ (A.6) 

where Ad^B = A(d fl B) — (d fl A)B. Finally, the FFS (fermion-fermion-scalar) 
couplings, ga lFiS \ are defined by the following Lagrangian term: 

C = g^ s 7\P a F 2 S, (A.7) 

where a = L or R. These simple rules determine uniquely the magnitude and the 
phases of the MSSM couplings that appear in Appendices and |D[ 

B Two point functions 

B.l Oblique parameters S,T,U,R 

The transverse parts of the four gauge-boson-propagator functions are parametrized 



in ref. 123 



as 



rr?V) 


= e n T (q ), 




(B.la) 


nfV) 


= e#z{r4 Q (g 2 ) 




(B.lb) 


nf(g 2 ) 


= 


-25 2 n?(g 2 )+5 4 n? Q (g 2 )}, 


(B.lc) 


IT T (g 2 ) 






(B.ld) 
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where g 2 (/x) = e 2 (/i) / s 2 (/i) = g 2 z {ji)c 2 {jj) are the MS couplings and all the polar- 
ization functions are renormalized in the MS scheme. The four effective charges 
are defined by 



1 


i 


e 2 (q 2 ) 


e 2 (/i) 


*V) 




e 2 (q 2 ) 


e 2 (/i) 


1 


f 






1 


1 







+ ReT^% 2 



ReU 3 T Q (q 2 ), 



+ Ren 



- 2s 2 ReIT^(g 2 ) + s^ReU^ z {q 2 



T,Z\ 



+ Ren T w (q 



where 



■=AB 



n 



t,v(q ) 



U-(q 2 



n 



T 



[m 



v) 



q 2 — my 



(B.2a) 
(B.2b) 
(B.2c) 
(B.2d) 

(B.3) 



and 'overline' denotes the inclusion of the pinch term g2| . The oblique correction 
terms are then expressed in terms of these effective charges and the weak boson 
masses 



23 



S 
4 

aT 

S + U 
4 

R 



s 2 (m 2 z )c 2 (m 2 z ) 



a{m z 



Aix 

qWY 



m 2 w #f(0) : 
s 2 (m|) 47r 



167T 



9 2 W W 

167T 



$1(0) 



or of the reduced functions, TT^ v (g 2 ), IT^ v (g 2 ), irj?(g 2 ) and IT^ 1 (q 2 ) 

w33 



5 

T 
U 
R 



167rRe 
4y/2G F 



n;%(o) 



16vrRe 



n?(o) -n"(o) 



4 3 z(o) 



167T 



n?>) - < z (g 2 ) - 2S 2 {n^(0) - if T Q z ( q 2 
+ 5 4 {n?S(o)-ng(g 2 )}. 



(B.4a) 
(B.4b) 
(B.4c) 
(B.4d) 

(B.5a) 

(B.5b) 
(B.5c) 

(B.5d) 



The MSSM contribution to the S, T, U, R parameters are then calculable from the 
MSSM particle contributions to the four IT^ 5 functions. All the results listed 



below agree with those in ref. |26 
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B.2 MSSM Higgs bosons 

The MSSM Higgs boson contributions are summarized as follows: 
167r 2 n« Q = B 5 (q 2 :m H -,m H -), 

167r 2 n^ Q = 7}B 5 (q 2 : m H -,m H -), 



(B.6a) 
(B.6b) 



i67r 2 n3? 



cos 2 (ci! - [3){B 5 (q 2 : m h ,m A ) + B 5 (q 2 : m z ,m H )} 



+ sin 2 (a - (3){B b (q 2 : m H ,m A ) + B 5 (q 2 : m z ,m h )} + B 5 (q 2 : m H -,m H - 



+ ™ 2 z 



i67rn 



2tt11 



sin 2 (a — P)B (q 2 : m z , rrih) + cos 2 (a — f3)B (q 2 : m z , niu) 
cos 2 (a - (3){B b (q 2 : m H -,m h ) + B b (q 2 : m w , m H )} 



(B.6c) 



+ sin 2 (a - (3){B b (q 2 : m H -,m H ) + B b (q 2 : m w ,m h )} + B 5 (q 2 : m A ,m H - 



+ m 



IV 



sin 2 (a - P)B Q (q : m ft ) + cos 2 (a - f3)B (q : m w , m fl ) 



(B.6d) 



The i?-functions are defined in Appendix D in ref. p3| . The SM Higgs boson 
contribution is reproduced by removing all terms with H~,H and A, and by 
setting sin 2 (a — (3) — 1 and rrih = rn>H SM - 

B.3 Scalar fermions 

The scalar fermion contributions are summarized as: 

2 

r 2rrQQ 

a=l 

2 

r 2rr3<2 

a=l 

2 



i67r 2 n^ g = c>qJ J2 B ^i 2 '■ m f a i m fji 

a=l 

167r 2 n^ = C f Q f h f £ \(Uf) al \ 2 B 5 (q 2 : mj^mjj, 

a=l 

16n 2 Uf = C f I 2 f £ |(f/ 7 )airi(t/V| 2 i?5(g 2 :m 7Q ,m 7 ; 



(B.7a) 
(B.7b) 
(B.7c) 



a,/3=l 



16vr 2 H^ = C,- £ K^iH^Vl W : m~ a ,m^), (B.7d) 

a,/3=l 

where the color factor Cf is 3 for squarks and 1 for sleptons. The electric charge 
Qf is given by (2/3, —1/3,0, —1) for (u,d,Pi,l). I%f denotes the third component 
of the weak isospin: +1/2, —1/2 for the up- and down-type sfermions, respectively. 
The indices a, (3 take L,R for the first two generations (no left-right mixing) and 
1, 2 for the third generation. If there is no mixing between fi and fn, the unitary 
matrix (U*) a p should be replaced by 5 a p. 
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B.4 Charginos and Neutralinos 

The chargino and neutralino contributions are as follows: 

16vr 2 n^ Q = 8q 2 B 3 (q 2 : m~-, m~ 

2- 



2 



^2q 2 B 3 (q 2 : m~~,m~ 



(B.8a) 
(B.8b) 



167r 2 n^ 3 



2{|pL) y f + |pji) -| 2 }(2 g 2 S3 ~ 5 4 ) 



(g 2 :m~ m~ ) 



+ 



^N^ 2 + \{N R ) l3 \ 2 y2q 2 B, - B A ) 



+m~ m~ {Nl^Nr)^ + (N L )*AN R )ij \B 



(q 2 : m~o,m~o),(B.8c) 



i67rn 



2tt11 



|(Cx)<il + 1(^)^1 \(WB 3 - B 4 )(q 2 : m~o,m~-) 



+m^m~-\{C L ) ij {C R )* + {Ci)UCr)h \B (q 2 : m~ ,m~-) 



where 



(B.8d) 

(B.9a) 
(B.9b) 
(B.9c) 



Here a = L or i? in (|B.9a|) and (|B.9b|) 
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C Vertex corrections on Z 



fa fa 



Here we give the 1-loop corrections to the process Z — > f a f a in the MSSM. The 
radiative corrections to the effective coupling g£ is denoted by 



*9i 



4V2G F m 2 z 



^ /z s' /Q (o)-r /Qi 



rrir 



(CI) 



where E^ a (0) is the derivative of the self energy function of the external fermion 
f a with the chirality a = L or R, whose mass is neglected: 



d 



*/.«» = 



q 2 =0 



(C.2) 



We can express the self energy E^(0) and the vertex function Tf a (q 2 ) mediated 
by a fermion ip and a scalar in a compact generic notation as follows: 



(4tt) 2 s' /q (o) = q 

(4vr) 2 r /tt (g 2 ) = 



B Q + Bi (0 : m^,m^), 



(C.3a) 



-g 2 (Ci2 + C 23 )-2C 2 4 + - 



x (p!,p2 : m^,m 0fc ,m^.) 



.(C.3b) 



Here C 9 is 4/3 for the gluino contribution (■?/> = g) and 1 for the others. The 
chirality index —a follows the rule: — L = R, —R = L. For each external fermion 
/, the following combination of {ip, 0} contribute to the vertices: 









/ = « 


d 




I 




= chargino 


{x + ,d*} 


{x;u*} 


{x + jn 




V> 


= neutralino 


{X°u*} 


{X°d*} 


{x°m 


{x%} (c 


i> 


= gluino 


{9,u*} 


{g,d*} 






4> 


= charged Higgs 


{d,H-} 


{u,H+} 


{l,H~} 


{u,H+} 





= neutral Higgs 


{it, (h, H, A)} 


{d,(h,H,A)} 




{l,(h,H,A)} 



The generic coupling notations of ( |A.5| ) - (|A.7|) then suffice to calculate all the 
Zfafa vertex corrections, Ag^. Summation should be taken over all non- vanishing 
coupling combinations. When d = b or i = r, the summations over the sfermion 
index is taken over i — 1 and 2. Otherwise, only the i = a sfermion contributes. 
We reproduced numerically the results reported in refs. [§, 0, |Kj . 
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D SUSY contributions to the /i-decay process 

The //-decay constant is parametrized in ref. [2^] as 

G p 



f w (0)+g 2 5 G 



AV2m 2 



w 



with 



D.l) 



(D.2) 



In the MSSM, A5 G receives contributions from the vertex and the box diagrams 



A6 G = 26^ + 5i b \ 



(D.3) 



D.l Vertex corrections 

The /1/2W vertex correction can be expressed as 

g hW 5 t«) = r W (0) _ ^^{^(o) + S' /2 (0)}. 

In the \i — > z/^ez/g process, the fiu^W' and v e eW + vertices give the some 5^ v \ with 
= g v ^ eW+ = JL. From the u e eW + vertex, we find 



:D.r 



9l 



Xi^e^e 

9l 



( J B + 5 1 )(0:m~ i ,m~o) + 



(B Q + B 1 )(0:m~ e ,m~o) + 



9l 



9l 



(fl o + £i)(0:m~,m~ 7 ), 
(D.2a) 

(D.2b) 



(^47Tj 1 eu e W+ \9l ) 9l \9l m x° m xr^o+ 9r 



X ^ (-2c 24 + -: 



x (0:m^,m~ ,m~-) 



\9l ) 9l \9l m x° i m xj Co + 9r " 



Y?Y~W 1 

-2C 24 + - 



x (0 : m~ r ,m~ e ,m~o) 



+ lfi , L J 9l ^24(U . m~ ,m~o,m~). 



In eqs. (p.2| ), ( p.3| ), summation over i = 1 to 4 (%°) and j = 1 to 2 (%j) should 
be understood. The momentum arguments of the C-functions are set to zero 

(p? = p*Pj = 0). 
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D.2 Box corrections 

The box contributions to the jj, — > v^eve amplitude can be expressed as 

iT = z{m(1) + M(2) + M(3) + M^W^PlVvJU^^Plu^ (D.4) 

By taking into account the normalization of the tree-level amplitude, — ^ 2 /2m^, 
the box diagram contributions to the 5 G parameter is 

i^-^VMfi). (D.5) 
9 2 £1 

Each M(i) is given by 
16tt M(l) = )^ (# L J ) D 2r (m~ L ,m 7L ,m~+,m~o) 

(D.6a) 

16tt M (2) = ) c/ L J (# L l ) p L « D 27 (m~ /i ,m~,m~-,TO~o) 

(D.6b) 

16vr 2 M(3) = ^jP^^)1^)'5oK^,^, m? ) 

(D.6c) 

16tt M(4) = -m~om~-g L ' g L l (g L 3 ) (g L J ) D (m~ , m~ m~-, m~o) 



2 x » x 



(D.6d) 

All the D-functions are evaluated at the zero momentum transfer limit. We re- 



produced the results presented in ref. [43 . 
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